Introduction to Nonequilibrium Statistical Mechanics 
with Quantum Field Theory 

Takafumi KiTA 
Department of Physics, Hokkaido University, Sapporo 060-0810, Japan 

^ . (Received December 30, 2009; revised February 27, 2010) 

^ ' In this article, we present a concise and self-contained introduction to nonequilibrium 

^-^, statistical mechanics with quantum field theory by considering an ensemble of interacting 

C^ . identical bosons or fermions as an example. Readers are assumed to be familiar with the 

Matsubara formalism of equilibrium statistical mechanics such as Feynman diagrams, the 
proper self-energy, and Dyson's equation. The aims are threefold: (i) to explain the fun- 

^f) , damentals of nonequilibrium quantum field theory as simple as possible on the basis of 

the knowledge of the equilibrium counterpart; (ii) to elucidate the hierarchy in describing 
nonequilibrium systems from Dyson's equation on the Keldysh contour to the Navier-Stokes 

1-^ I equation in fluid mechanics via quantum transport equations and the Boltzmann equation; 

(iii) to derive an expression of nonequilibrium entropy that evolves with time. In stage (i), 
we introduce nonequilibrium Green's function and the self-energy uniquely on the round- 
trip Keldysh contour, thereby avoiding possible confusions that may arise from defining 
multiple Green's functions at the very beginning. We try to present the Feynman rules 
'j^ . for the perturbation expansion as simple as possible. In particular, we focus on the self- 

consistent perturbation expansion with the Luttinger-Ward thermodynamic functional, i.e., 
Baym's ^-derivable approximation, which has a crucial property for nonequilibrium sys- 
tems of obeying various conservation laws automatically. We also show how the two-particle 
correlations can be calculated within the ^-derivable approximation, i.e., an issue of how to 
handle the "Bogoliubov-Born-Green-Kirkwood-Yvons (BBGKY) hierarchy" . Aim (ii) is per- 
formed through successive reductions of relevant variables with the Wigner transformation, 
the gradient expansion based on the Groenewold-Moyal product, and Enskog's expansion 
from local equilibrium. This part may be helpful for convincing readers that nonequilibrium 
systems can be handled microscopically with quantum field theory, including fluctuations. 
We also discuss a derivation of the quantum transport equations for electrons in electromag- 
netic flelds based on the gauge-invariant Wigner transformation so that the Lorentz force 
is reproduced naturally. As for (iii), the Gibbs entropy of equilibrium statistical mechanics 
^ ' suffers from the flaw that it does not evolve in time. We show here that a microscopic ex- 

CO , pression of nonequilibrium dynamical entropy can be derived from the quantum transport 

On ■ equations so as to be compatible with the law of increase in entropy as well as equilibrium 

CO ' statistical mechanics. 
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§1. Introduction 

Nonequilibrium quantum field theory on the real-time Keldysh contoupl^ has 
been used extensively in recent years to describe a wide range of dynamical phe- 
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nomena in condensed matter physics, nuclear matter, and high-energy physics.'^''^ 
The method enables us to handle a wide range of many-body systems with quan- 
tum effects and/or strong correlations microscopically from first principles. There 
are already excellent review articlea^^^ and textbooks^'l^ on the topic. The ap- 
proach is also explained briefly in some standard textbooks on statistical many-body 
theory.l^'l^ Thus, adding another review here on the topic may require some jus- 
tifications for doing so. The distinct features of this article may be summarized as 
follows. 

(i) Considering specifically a system of identical bosons or fermions with a two- 
body interaction, we try to explain the basics of nonequilibrium quantum field theory 
concisely and clearly on the basis of the knowledge of the equilibrium counterpart. 
For example. Green's function and the self-energy are introduced uniquely on the 
round-trip Keldysh contour in exactly the same way as those on the imaginary-time 
Matsubara contour. We thereby avoid defining multiple real-time Green's functions 
at the very beginning so as not to cause confusion. We also try to present Feynman 
rules for the perturbation expansion as simple as possible. One of the extra factors 
here is how to carry out the summations over the forward and backward paths of 
the Keldysh contour. Several variant rules have been proposed for it such as the one 
by KeldyshP and another by Langreth (i.e., the Langreth theorem).'^ We present 
here an alternative rule for the summation, which we think is the simplest possible. 
It is worth pointing out that the consideration here can be extended easily to other 
systems, e.g., electrons in solids with a periodic potential, impurities, phonons, etc. 
Indeed, we only need to modify the basic Hamiltonian to incorporate those effects. 
This may be realized by looking at the variety of applications of the method ranging 
from condensed matter physics to high-energy physics .1^^^ 

(ii) We specifically focus on the self-consistent perturbation expansion with 
the Luttinger-Ward thermodynamic functional,'^ which reproduces the Hartree- 
Fock theory as the lowest-order approximation.l^ Its essence lies in determining 
Green's function G and the self-energy E self-consistently using Dyson's equation 
{Gq — Z!)G = 1 and the relation S oc 6^/6G based on a given functional <P = ^[G]. 
As shown by Baym,!^ the approximation scheme has a unique property of obeying 
various conservation laws automatically, i.e., a property indispensable for describ- 
ing nonequilibrium systems but not satisfied by the simple perturbation expansion. 
Despite this crucial feature, the self-consistent perturbation expansion has not been 
paid sufficient attention in those review articles and textbooks. We give it a full 
account together with a complete proof of the conservation laws. We also derive 
the Bethe-Salpeter equation for two-particle correlation functions in the (^-derivable 
approximation; thus, choosing a definite ^ = ^[G] will be shown to amount to deter- 
mining the whole Bogoliubov-Born-Green-Kirkwood-Yvons (BBGKY) hierarchy.^^ 

(iii) We try to elucidate the hierarchy in describing nonequilibrium systems from 
the microscopic Dyson's equation on the Keldysh contour to the macroscopic Navier- 
Stokes equation; the latter is shown to be reached from the former by successive 
reductions of relevant variables through quantum transport equations and the Boltz- 
mann equation. The Navier-Stokes equation is an archetype of nonlinear evolution 
equations on which a wide variety of nonequilibrium phenomena (e.g., turbulence. 
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chaos, and pattern formation) have been discussed extensively!^''^ A manifest 
derivation of it from Dyson's equation may convince readers that those phenomena 
can be handled microscopically from first principles and enable them to incorporate 
effects beyond the phenomenological approach such as fluctuations. 

(iv) We derive a microscopic expression of nonequilibrium entropy, which evolves 
with time. Entropy S is the central quantity in thermodynamics as embodied in the 
Clausius inequality dS > TdQ, where T and Q denote temperature and heat, respec- 
tively. It tells us that entropy for any isolated system should increase monotonically 
with time. On the other hand, it was shown by Jaynea^ that every equilibrium 
statistical ensemble can be identified as the maximum of the Gibbs (or information) 
entropy 

Y^ , I ks '■ Boltzmann constant 



S- kB^^Pulnp,, I : probability of the state z. ' ^^'^^ 

under some constraints. However, the Gibbs entropy suffers from the lack of dynam- 
ics, i.e., it cannot describe the law of increase in entropy.'^ Except for the Boltz- 
mann entropy for dilute classical gases^ '' ' and its extensions to classical denser 
systems,'^''^ we do not have a widely accepted statistical-mechanical expression 
for nonequilibrium entropy to represent the Clausius inequality microscopically. Fol- 
lowing Ivanov et al.,'^ we show here that such an expression can be derived from 
the quantum transport equations mentioned above so as to be compatible with the 
law of increase in entropy and also to embrace the Boltzmann entropy as a limit. 

(v) We extend the derivation of quantum transport equations to electrons in 
electromagnetic fields, where special care is necessary for the gauge invariance of the 
equation to reproduce the Lorentz force adequately. 

(vi) To supplement contents above, we include in Appendix A a full account of 
the second quantization method as an equivalent alternative to the description with 
many-body wave functions in the configuration space. We also describe the equilib- 
rium self-consistent perturbation expansion with ^ in Appendix B, the Luttinger- 
Ward thermodynamic functional in Appendix C, and a derivation of an expression 
of equilibrium entropy in Appendix D. 

Now, we briefly summarize the relevant history of quantum field theory in statis- 
tical mechanics. The application of the method to equilibrium many-body phenom- 
ena was pioneered by Matsubara in 1955 based on the simple perturbation expan- 
sion with the imaginary-time Matsubara Green's function. ■^^^ The finite-temperature 
version of Wick's theorem,-—' which forms the basis for the perturbation expansion 
with the Matsubara Green's function, was first proved by Bloch and de Domini- 
cia^ in 1959 and refined by GaudinPSJ shortly. In 1960, Luttinger and WarcP^ 
provided a formally exact expression of the thermodynamic potential as a func- 
tional of the fully renormalized Green's function G, which contains the functional 
^[G] mentioned above in the form of the skeleton diagram expansion. It was subse- 
quently used by Luttingepl^ to clarify some general properties of interacting normal 
fermions, e.g., the Fermi-surface sum rule anticipated by Landau in his Fermi-liquid 
theory.l^ Luttinger^ also noted that the skeleton diagram expansion forms a ba- 
sis for a self-consistent perturbation expansion where the Hartree-Fock theory is 
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reproduced as the lowest-order approximation. Shortly, equilibrium quantum field 
theory was made quite popular among theorists through the well-known textbook 
by Abrikosov, Gor'kov, and DzyaloshinskipO which even contains its applications 
to superconductivity and Bose-Einstein condensation as well as a microscopic justi- 
fication of the Landau Fermi-liquid theory. Meanwhile, the method had been used 
fairly extensively for describing nonequilibrium transport phenomena by tracing the 
time evolution of the density matrix.'^^ A crucial step was made along this line 
by KadanofF and Baym in 1962,'^ who pioneered the use of Green's function for 
this purpose. However, their treatment was still based on the analytic continuation 
from the imaginary-time Matsubara Green's function. BaynpSl' successively clari- 
fied a sufficient condition for various conservation laws to be satisfied automatically 
in terms of 'P[G] and called it "^-derivable approximation". The approximation 
scheme is essentially identical in content with the self-consistent approximation of 
Luttinger but superior to the latter in its explicit usage of ^\^ On the other hand, 
Schwinger^ introduced Green's functions on a closed time-path contour in 1961 to 
calculate time-dependent expectation values of physical quantities for a Brownian 
motion. In 1964, KeldysbpJ'''^ clearly recognized that quantum field theory with 
Green's function can treat general nonequilibrium phenomena in the same way as 
equilibrium phenomena to develop a perturbation expansion with respect to the in- 
teraction on the real-time round-trip contour of — oo <t< oo. The theory made an 
epoch in nonequilibrium statistical mechanics with quantum field theory. However, 
its great potential for describing nonequilibrium phenomena seems to have been ap- 
preciated rather slowly. Some of its earliest applications include: a calculation of 
the tunneling current through a metal-insulator-metal junction by Caroli et al.j^ a 
derivation of a transport equation for a superconductor by Aronov and Grevichp^ a 
derivation of nonequilibrium quasiclassical equations of superconductivity by Larkin 
and Ovchinnikov.^^' Later developments may be found in Refs. [^- Ill|) . 

§2. Representations and expectation value 

We first introduce the Schrodinger, Heisenberg, and interaction representations 
for a time-dependent Hamiltonian. We then derive an expression for expectation 
values of physical quantities in terms of the latter two representations. 

2.1. Schrodinger and Heisenberg representations 

The Schrodinger and Heisenberg representations are usually introduced in terms 
of a time-independent Hamiltonian. Since we treat nonequilibrium dynamical phe- 
nomena, we start by defining those representations explicitly for a time-dependent 
Hamiltonian. 

We consider a many-body system described with a time-dependent Hamiltonian 
7i{t). By adopting the second quantization metnod, the Schrodinger equation reads 

^^^^^ = m\Mt)), (2-1) 

where \^i,{t)) denotes a many-body wave function specified by quantum number i/ 
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with (y/^(f)|!P'i,(t)) = 1. See Appendix A for the notations of second quantization as 
weh as for a proof that the method is completely equivalent to the description with 
many-body wave functions in the configuration space. Equation (|2-1|) has the formal 
solution: 

\Mt))=^it,to)\Mto)), (2-2) 



where Iq is some initial time and U{t,tQ) is defined by 



Uit,to) = l + J2 



n=l 



dtr. 



to 



t3 



t2 



dt2 / dtin{tn)---n{t2)n{ti). 

to J to 



(2-3a) 



It is seen easily that Eq. ([2^ with Eq. ([Ma]) satisfies Eq. ([20]) . Although t > to 
may have been assumed implicitly, definition ()2-3p is also effective for t < to- We 
then observe that, for t > to [t < to), the Hamiltonians on the right-hand side are 
arranged in chronological order from right to left (left to right). Hence, it follows 
that the operator U(t,to) can also be put into an exponential form as 



U{t,to) 



Texp 
Texp 



-- / n{ti)dti 

"-J to 



n{ti)dti 



to 



t>to 
t<to 



(2-3b) 



where T (T) is the time-ordering (anti-time-ordering) operator to arrange the Hamil- 
tonians into the chronological order from right to left (left to right) with an extra 
sign change upon every transposition of fermion field operatorsp21i''21i''l2SI' this sign 
change is irrelevant here as ?^(t) is usually composed of an even number of field 
operators. See, e.g., p. 57 of Ref. [39l) for the equivalence between Eq. (|2-3ap and Eq. 
(|2^ . 

One can show easily that U{t,to) obeys the equation of motion: 



ih 



dt 



n{t)u{tM)- 



(2-4a) 



Expressing the integrals in Eq. ()2-3p as /^ dti j^ dt2 • • • /j dtn, we also obtain 



ih 



d K{t,to) 
dtn 



-U{t,to)'Hito 



(2.4b) 



The operator has the properties: 

U{t,ti)Vl{ti,to)=llit,to), (2-5a) 

Uito,to) = l. (2-5b) 

Equation ()2-5bp is obvious from Eq. ()2-3ap . whereas Eq. ()2-5ap can be proved as 
follows. We first observe that V(t,to) = ^{t,ti)IJ{ti,to) satisfies the same first- 
order differential equation with respect to t as IJ{t,to). We also notice with Eq. 
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Fig. 1. Integration contour of Eq. H2-7|) . Both the forward and backward paths are actually on the 
real axis but shifted slightly upwards and downwards, respectively, to distinguish them clearly. 

()2-5bp that the initial values at t = ti are the same between the two operators, i.e., 
V(t,to) =U{t,tQ). Hence, we arrive at Eq. (j2-5ap . 

Setting i = to in Eq. (|2-5a|) enables us to identify lA~^{t, to) = ^(to, t). Moreover, 
taking the Hermitian conjugate of Eq. ()2-4ap and comparing the result with Eq. 
()2-4bp . we obtain Z^^(t,to) = Z^(to,t). These results are summarized as 

U-\t,tQ)=U\tM)=U{to,t). (2-6) 

With lA, the expectation value of an Hermitian operator 0{t) in terms of wave 
function (|2-2p can be expressed alternatively as 

{i^,{t)\d{t)\K{t)) = {Mh)\dn{t)\Mto)), (2-7) 

where 0}i{t) denotes the Heisenberg representation: 

dn{t)=U\tM)d{t)U{tM)- (2-8) 

Note that C'H(t) is Hermitian with to dependence through U\ the latter dependence 
is dropped here for convenience. The left (right) side of Eq. (j2-7p corresponds to 
the Schrodinger (Heisenberg) picture. We realize with ^^(t, to) = U(to,t) that the 
evaluation of Eq. ()2-7p in the Heisenberg picture proceeds along the closed time path 
of Fig. 1. 

2.2. Interaction representation 

We now consider those cases where Hit) can be split into two parts as Ii{t) = 
Ii{){t) + T-L'{t). Following Eq. ()2-3p . we first define a time-evolution operator with 
Mt) by 

Uoit,to) = l + y^(--) dtn--- dt2 dtino{tn)---noit2)noiti), (2-9) 

which clearly has the same properties in terms of T-Lq as Eqs. (|2-4p - (|2-6p of li with 
%. We next introduce the so-called scattering matrix: 

5(t,to)=Z^o(to,t)^/(t,to). (2-10) 

Its inverse and Hermitian conjugate are easily identified with Eq. (|2-6p as 

S-\t, to) = S\t, to) = Zi(to, t)Zio(i, to)- (2-11) 

Note that 5~-'^(t,to) / S(to,t) generally. Using Eqs. (|2-4p and (|2-5p . one can show 
easily that 5(t,to) satisfies 

^h^^^=n[it)S{t,to), (2-12) 
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where 7^j(t) denotes the interaction representation of T-L'{t); the interaction repre- 
sentation is defined for an arbitrary Herniitian operator 0{t) by 

di{t) =Uo{tQ,t)d{t)Uo{t,to)- (2-13) 

Note that C'i(t) is also Hermitian with to dependence through Uq; the latter depen- 
dence is also dropped here for convenience. Noting Eqs. (|2-3ap and (j2-4ap . we can 
immediately write down the solution to Eq. (|2-12|) as 

S{t,to) = l + y^{--] dtn--- dt2 dh-H[{tn)---n[{t2)n'i{tl). (2-14) 

„_i V "-/ Jto Jto Jto 



It may also be put into an exponential form as Eq. ()2-3bp . 

Let us substitute Z^(t,to) = l(oit,to)S{t,to) from Eq. (120011 into Eq. ([23]). We 
thereby obtain an alternative expression for the expectation value in terms of Eqs. 
(|2T3]1 and ([2011) as 

{Mt)\OimAt)) = {Mto)\SHt,to)di{t)S{t,to)\^,ito)). (2-15) 

2.3. Keldysh contour 

From now on, we set tQ = —oo. Using S'^{oo,t)S{oo,t) = 1 and 5(t2,ii)5(ti, to) = 
5(^2,^0)) we can transform Eq. ()2-15p as 

{K{t)\d\K{t)) = (i^.(-oo)|5t(t,-oo)Oi(t)5(t,-oo)|'2^,(-oo)) 

= {^u{-oo)\S\oo,-oo)S{oo,t)Oi{t)S{t,-oo)\^u{-oo)) 

= {^y{-oo)\S\t,-oo)di{t)S\oo,t)S{oo,-oo)\^u{-oo)). (2-16) 



We also notice that Eq. (|2-14p with to = —00 is equivalent to 



5(t, —00) = Texp 



- / n[{t{)dh 



(2-17a) 



The Hermitian conjugate of Eq. ()2-14p can be expressed with the anti-time-ordering 
operator T as 

i 



5' (t, —00) = Texp 



'H[{ti)dti 
" -It 



(2-17b) 



Substituting these expressions into it, we realize that the evaluation of Eq. (j2-16p 
proceeds on the round-trip Keldysh contour C of Fig. 2 starting from t^ = —00 
towards t = 00, and C'i(t) can be evaluated on either of the forward or backward 
path equivalently. 

C 



t 
Fig. 2. Keldysh contour. 
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The fact just mentioned enables us to provide a concise expression to Eq. ()2-16p . 
Let us introduce the operator Sc on C by 



>c 



Tcexp 



■ym-)<it< 



(2-18) 



where t denotes a point on C, and Tc arranges operators according to their chrono- 
logical order on C from right to left. By using Sc, Eq. (j2-16p can be expressed as 

{•PAt)\6\m)) = {M-oo)\TcScd,mA-o.)) = (^-(-°°)l7-cfe6,W|i>.( oo)) ^ 

{'F^{-oo)\Sc\M-oo)) 

(2-19) 
where C'i(t) can be on either the forward or backward path, and we have used 
{'l^u{—(X))\Sc\'I^u{—co)) = 1 in the last equality. One can see clearly from the last 
expression that only those diagrams that are connected with Oi contribute to the 
expectation value. 

§3. Nonequilibrium perturbation expansion 

With the above preliminaries, we now embark on formulating nonequilibrium 
statistical mechanics for an ensemble of identical bosons or fermions with quantum 
field theory. There have been basically two approaches regarding the initial condi- 
tions at to = — oo. The first one is that of Keldysh,^- where the initial state is void of 
interaction to switch it on adiabatically together with external fields. The second one 
is the grand canonical ensemble with full interaction; however, the initial correlations 
along the imaginary-time path are neglected subsequently.'^''^ We adopt here the 
approach by Keldysh as it is free from any approximations and, hence, transparent. 
Indeed, the effect of the initial correlations may be handled using Keldysh's approach 
by waiting thermalization due to interaction before applying external fields. 

3.1. Hamiltonian, density matrix, and expectation value 

We consider a system of identical noninteracting bosons or fermions at to = — cx). 
Let us apply an external field U{r,t) for t > — oo and also introduce the interaction 
'Hint adiabatically. Thus, the total Hamiltonian at time t is given by 

■H{t)=-Ho{t)+Uir.ta{t). (3-1) 

Here, I-Lq denotes noninteracting Hamiltonian including U , limt describes the two- 
body interaction, and a{t) is some adiabatic factor. The last quantity may be ex- 
pressed with the step function 9{t) as a{t) = 9{—t)e^+^ + 9{t), for example, with 
0+ denoting an infinitesimal positive constant. By neglecting the spin degrees of 
freedom for a while, T^o and liuit can be expressed as 

7io= U\ri)ki{l,{ri)d\i, (3-2a) 



n^nt = \ I d\i I d\[^^{ri)^\r[)V{ri-r'Mr'Mri). (3-2b) 
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Here, Ki denotes 

Ki = -—V? + [/(nil) - /i, (3-3) 

Im 

with m and // as the particle mass and chemical potential, respectively; the field 

operators ijj and ■0' satisfy 

i){ri)'4}\r2) =F ■0^(^-2 )V'(n) = <J(n - 7-2), (3-4a) 

V'(ri)Vi(r2) =F V^(r-2)V'(n) = 0, (3-4b) 

with the upper (lower) sign corresponding to bosons (fermions), and V{ri — r^) is 
the interaction potential with V{—r) = V{r), which may be expanded as 

with V-p = Vp. 

Let us assume that the system at tg = "Oo is in thermodynamic equilibrium 
described with the grand canonical ensemble. The corresponding density matrix 
po can be expressed with the eigenvalue Ey{—<x)) and the eigenstate \^y{—<x))) of 
?^o(— oo) as 

Zq — ' — ' 

(3-6) 
where Tq and /cb are the initial temperature and Boltzmann constant, respectively. 
We now average Eq. ()2-7p over the initial probability distribution Z^ g^~'Eu{-oo)/k^TQ ^ 
Then, the expectation value of 0]i{t) can be expressed as 

(OhIO) - TrpoOH(t) = iZ^Z^I^, (3.7) 

where we have used Eq. (j2-19p in the second equality. Since the initial state is a 
noninteracting grand canonical ensemble, we can use the Bloch-de Dominicis theo- 
rem,'23'l22li i.e., the finite-temperature version of Wick's theorem,'^ to evaluate the 
second expression of Eq. p-7p perturbatively in terms of T^int- It then follows that 
only those Feynman diagrams connected with Oi{t) contribute to the expectation 
value due to the denominator in Eq. (|3-7p . 

A few comments are in order. First, the original postulate of the noninteracting 
grand canonical ensemble in proving the Bloch-de Dominicis theorem can be relaxed, 
as already pointed out by Danielewicz.'SJ To be specific, the Wick decomposition is 
justified for any density matrix of the form: 

'5o = nEl^'=)Pfc'^"'^l' l^fc) = H=|0)a:, (3-8) 



where pk denotes the probability that there is a particle in the single-particle state 
k, Uf^ is the number of particles occupying k, c|. is the creation operator, and \0)k is 
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the vacuum of k; see also Appendix I A. 91 for the equihbrium density matrix. This po 
is easily shown to satisfy 

ft 
CkPo=PkPoCk, Cf^Po = — PoCk- (3-9) 

Pk 

Using these relations, one may repeat the concise proof by Gaudiii^''^ to convince 
oneself that the Bloch-de Dominicis theorem holds even for pk ^ Q-i'^k-N/kBTo^ where 
ek is the single-particle energy of the state k. Thus, we can evolve the system also 
from noninter acting nonequilibrium states of a wide class. 

Second, the effects of so-called "initial correlations" may be incorporated into 
the formalism by waiting for thermalization due to interaction before applying exter- 
nal fields. Indeed, thermalization is achieved in this formalism by considering terms 
of the second (and higher) order in the perturbation expansion, which describe colli- 
sions between particles. In this context, it may seem puzzling why thermalization is 
achieved only by evolving the system mechanically. Equivalently, what is the origin 
of the arrow of time in this formalism? It is instructive for answering the question to 
see the results of a molecular dynamics simulation on a classical dilute hard-sphere 
gas by Orban and Bellemans.^'^' They show that the entropy of the system (i.e., 
the Boltzmann entropy that is applicable here) does increase for the overwhelming 
majority of initial conditions without any contradictions to the time-reversal sym- 
metry in the mechanical laws. Suggested by the results and following the viewpoint 
of Boltzmann,-^ we trace the origin of the arrow of time in many-body systems to 
the overwhelmingness of those initial conditions. The equation of motion for the 
nonequilibrium Green's function below will be regarded as corresponding to those 
overwhelming majority of initial conditions. Thermalization will also be achieved 
for systems in contact with thermal baths. The concept of temperature is not nec- 
essary at all in this formalism; it naturally emerges as one of the properties in the 
thermalized distribution function. 

3.2. Green's function 

We now consider the effects of the interaction 

n'it) = Hintait) (3-10) 

perturbatively. A key quantity to this end is Green's function: 

G(l^,2^) = -l(rcViH(l^)V'^(2^)> = -'-{TcSc4'ii^'')4,li2^)),. (3-11) 

Here, the superscript '^ in 1*-^ = ritf distinguishes the forward and backward paths 
of the Keldysh contour, and the subscript c in the last expression denotes retaining 
only connected diagrams. We drop the subscript i in ■0i and t/^l below; we also omit 
the superscript *-^ whenever it is unnecessary such as those in the argument of V. 

The perturbation expansion on C can be carried out in the same manner as 
that of the equilibrium Matsubara formalism on the imaginary-time contour t S 
[0, —ih/k^T] with the difference lying only in the contours; see Appendix |B] for the 
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Matsubara formalism. Hence, it follows that Eq. ()3-lip satisfies Dyson's equation: 

^^] (13^ = 6(1^,2^), (3-12) 



(ih^ - k^^ G(l^, 2^) - I i:(l^, 3^)G(3^, 2^ 



where Ki is given by Eq. (j3-3p . 

Following Keldysh/' we next transform every integration on C into that of the 
ordinary time t £ [— oo, oo]. To this end, we express C as a sum of the forward path 
Ci and the backward path C2 as 



[ dt^ = I dt^ + f df = f dt^ - I 

Jc JCi JC2 i-oo(Ci) J- 



dt' 



-oo(C2) 



(3-13) 



We next introduce 1* = rit\ (i = 1, 2) to distinguish times on Ci and C2 clearly. We 
then define G*-'(l,2) = G(l*,2-') to provide them with the matrix representation: 

" G"(l,2) Gi2(l,2) " 
G21(1,2) G22(1,2) 



6(1,2)^ 



(3-14) 



We express every 2x2 matrix due to the path decomposition with " on top of it like 
Eq. ()3-14p . Explicit expressions of G*-' are given by 



G^^(l,2 

il2 



-^(V'h(1)^1i(2)), 



G-(l,2)=T^(ViJi(2)^H(l)), 



G''{1,2) 



G^\l,2) 



[9{h-t2){i^n{l)i'lm ± 9{t2-h){4;l{2)i;n{l))] 



h 

e{h-t2)G^Hl, 2) + 9{t2-h)G^\l, 2), 

-^[9{t2-ti){i,n{l)i'li{2))±9{ti-t2){4^l{2)i;n{l))] 
9{t2-ti)G^Hh2) + 9{ti-t2)G'\l,2), 



with 



9it) 



1 :t > 
:t < 



(3-15a) 
(3-15b) 

(3-15c) 

(3-15d) 
(3-16) 



The off-diagonal elements G^^ and G^^ correspond to G^ and G^ of Kadanoff and 
Baym,'^ respectively. The four elements have the properties: 



12/ 



^21/ 



[G^^(l,2)]* = -G^^(2,l), [G^^(l,2 



-<ii 



22/ 



G'^(2,l), [G^^(l,2)]* = -G'"(2,l), 



G^\l, 2) + G22(1, 2) = G^'^il, 2) + G2^(1, 2). 
Equation ()3-17p can be expressed concisely in terms of G as 

G(l,2) = -fiGt(2,l)fi, TYG(l,2)=TYfiG(l,2), 



(3-17a) 
(3-17b) 



(3-18) 



where fj (i = 1,2,3) are the Pauli matrices. It is also convenient to define the 
self-energy matrix: 

" i:ii(l,2) ^7^=^(1,2) " 

i;2i(i,2) i;22(i,2) 



r(i,2) = 



(3-19) 
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Now, we can express Dyson's equation ()3-12p in terms of G and IJ as 

ih^^ - k^ (5(1, 2) - / r(l, 3)f3G(3, 2) d3 = f35(l, 2). (3-20) 

Here, fs between G and 17 results from the transformation: 

/ U{l\3^)Gi3^,2^)d3^ = f d3[S'\l,3)G^^ (3,2) - S'^{1,3)G^^ {3,2)] , 



and fs on the right-hand side of Eq. (|3-20p originates from the opposite signs in the 
arguments of 9 in Eqs. (|3-15c|) and ()3-15dp . 
Let us introduce 



G(°)(l,2)=(.nJ--^i) f3(5(l,2), 



(3-21) 



which is the solution to Eq. ()3-20p without IJ, i.e.. Green's function without inter- 
action. Multiplying it by G^^'f^ from the right, Eq. ()3-20p can also be put into an 
integral form as 

G(1,2) = G(°)(1,2)+ 1^3 /"(i4G'W(l,3)f3^(3,4)f3G'(4,2). (3-22) 

3.3. Perturbation expansion of G 

Let us write down the Feynman rules for the perturbation expansion of G to 
see how they can be obtained from those of the equilibrium Matsubara formalism. 
Using Eq. (|3-13p , we first express Sc of Eq. (|2-18p as 



Sc = Tc exp 
= 7c'exp 



j=i 



^(-irW di^,n[{t 



- / dl / d2j;(-iriF(l-2)Vit(r)^t(2^)^(2^)^(r) 



i=l 



where V is defined by 



V{l-2) = 5{ti-t2)V{ri-r2). 



, (3-23) 



(3-24) 



This interaction can be expressed graphically as Fig. 3, where the wavy line denotes 
{—iy~^V and the two incoming (outgoing) arrows signify if) (V'^)- 




Fig. 3. Feynman diagram for the two-body interaction. 
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Let us compare Eq. (|3-23p with the corresponding S-matrix of the equihbrium 
Matsubara formaUsmpSJ 



5m(/3,0) ^T^exp 



I I .. I , .s rt^1^,/,t/ 



2^ , dl / d2y(l-2)^T(i)^T(2)^(2)^(i^ 



(3-25) 



where /3 = I/Zcb^ with T the temperature, every time integration extends over t G 
[0, —ih(3] with the periodic (antiperiodic) boundary condition for bosons (fermions), 
and Tt denotes the "time" -ordering operator on the imaginary contour; see also Eq. 
()B-3p of Appendix[B]for the derivation of Eq. ()3-25p . We then notice that, besides the 
difference in the integration contours, Eq. (j3-23p contains an additional summation 
over the paths i = 1, 2 with the potential (— 1)*~^V"(1 — 2). 

This observation enables us to obtain the Feynman rules for the nth-order con- 
tribution (n = 1, 2, • • • ) to G straightforwardly from those of the Matsubara formal- 
ism. 1S''3D'E3:IS1J They are summarized as follows: 

(a) Draw all possible nth-order connected diagrams. 

(b) With each such diagram, associate a factor 

^ ' ^ ^ (3-26) 



2"n! 

where n^ denotes the number of closed particle loops, and the factor (ih)'^ results 
from multiplying {—i/h)^ from the expansion of Eq. (|3-23p by (ih)'^^; the latter 
cancels each —i/h in definition ()3-lip of Green's function. 

(c) For an interaction line on the contour Cj (i = 1,2), associate the potential 

i-iy-^v. 

(d) With a particle line arriving at V from 2-', associate G*-''''''(l, 2). If the line 
enters and leaves the same interaction line with i = j, associate ^^^'''•'(l, 1+) 
and G^^^'^^(l+, 1) for i = 1,2, respectively, where the subscript + denotes the 
presence of an extra infinitesimal positive constant in the relevant time argu- 
ment. This rule works to recover the operator arrangement of Eq. (|3-2bp that 
V't lies to the left of V'- Note that G^^^'^^l, 1+) = G22(o)(1+, 1) = G^'^^^\l, 1), 
as seen from Eq. (j3-15p . 

(e) Integrate or sum over all the internal arguments. 

(f) Spin degrees of freedom may be incorporated by multiplying the contribution 
of each closed particle line by 25+1, where 5 denotes the magnitude of spin. 

It is worth pointing out that the factor (2"n!)~^ in Eq. ()3-26p is canceled eventually 
by the number of topologically identical diagrams in the perturbation expansion for 
G. 

Thus, one of the features inherent in the real-time nonequilibrium formalism lies 
in the additional summation over the paths Ci and C2. Several alternative versions 
of the Feynman rules have been presented for carrying out this extra summation. 
The Langreth theoreiEP'^^li ig among such alternatives. However, this rule with the 
path decomposition C — t- Ci + C2 may be the most natural, straightforward, and 
simple. 

As an example of these rules, let us consider the first-order contributions, where 
topologically distinct Feynman diagrams are given in Fig. 4. Those obtained with 



Introduction to Nonequilibrium. Statistical Mechanics 



15 




12 1' 12 



2' r 



Fig. 4. Topologically distinct first-order diagrams for G. Path indices are omitted. Those obtained 
from above with 2 •<->■ 2' yield the same contribution. 



the replacement 2 2' yield the same contribution to cancel the factor 1/(2^1!) 
in rule (b). Performing the summation over the paths Cj {i = 1,2), we obtain an 
analytic expression for the first-order Green's function as 

2 

xGi*'{0)(2, 1') + G*J'W(1,2)G^'^'W(2,2')G^'^'(°^(2', l')] 



^ / d2 / d2' G^^W (1, 2)5jy{-iy^^ih 
jj' 



±5(2,2') d3V{2-3) 



:G"^°\3,3) + y(2-2')G^'^'(°)(2,2') 



G^"*'(°)(2',l'), 



(3-27) 



where, in the second expression, we have factored out the two Green's functions 
connected with the external lines. A comparison of Eq. (|3-27p with Eq. (|3-22p leads 
to the expression of the first-order self-energy: 



^ij(si)(^^2) = 6iji-iy-Uh 



±(5(1, 2) Idl' y(l-l')G"(°) (1', 1') + y (1-2)G"(°) (1, 2 

(3-28a) 
where the superscript '^^ signifies that the expression belongs to the "simple" first- 
order perturbation expansion. Thus, the first-order self-energy is diagonal in the 
Keldysh space. We further consider the Feynman rule (d) above to make the replace- 
ments: G^W (i^ i') ^ g^W (i, i^) = ^i^W (i, i') and G^^W (1, 1') ^ G^^W (1+, 1') = 

G^'^^^'{1, 1'). We thereby obtain an expression for the first-order self-energy matrix 
r('^i) as 

^('^i'(l,2)=f3i;'^H^(l,2), (3-28b) 

with 
U'^^{1,2) = ±6{1,2) f dl' V{l-l')ihG^'^^°\l' ,1') + V{l-2)ihG^^'^'^\l,2). (3-29) 

We present further examples of the perturbation expansion with the Feynman rules 
in the next subsection. 

3.4. Self- consistent perturbation expansion for G 

The procedure described in §3.3 enables us to carry out the simple (bare) per- 
turbation expansion concisely up to a desired order or to sum up some partial series. 
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There is another category in the perturbation expansion, i.e., the self-consistent 
perturbation expansion, which is characterized by a unique property for describing 
nonequihbrium dynamical phenomena, i.e., various conservation laws are satisfied 
automatically. It was pioneered for a normal Fermi system by Luttinger and WavS^ 
in their effort to write down its thermodynamic potential rigorously as a functional 
of the exact self-energy U. Luttingerl^' subsequently pointed out that the expansion 
reproduces the Hartree-Fock approximation as the lowest-order approximation and 
may be used as a general approximation scheme. Baym™ ^ later presented a sufficient 
condition for various conservation laws to be obeyed automatically and showed that 
the self-consistent perturbation expansion based on the Luttinger- Ward functional 
meets the condition. The practical approximation scheme Baym presented there is 
the same as the one suggested by Luttinger^^ but superior to the latter in that it 
is given explicitly in terms of the functional ^. Indeed, it is with ^ that the conser- 
vation laws are proved; see §8 on this point. Moreover, the functional enables us to 
calculate two-particle and higher-order correlations systematically, as shown in §7. 

It is more convenient to express the Luttinger- Ward thermodynamic potential as 
a functional of Green's function GW^ Its key ingredient is the functional (p = (p[G] 
given as a power series in terms of the interaction potential, i.e., the so-called "skele- 
ton diagram expansion". Using <?, we can also construct self-consistent approxi- 
mations systematically up to a desired order. We now explain this approximation 
scheme in detail with its advantages together with further examples of the pertur- 
bation expansion with the Feynman rules. See Appendices [B] and O for the self- 
consistent perturbation expansion and the Luttinger- Ward functional, respectively, 
in the equilibrium Matsubara formalism. Although we limit our consideration below 
to normal phases, it is worth pointing out that the Luttinger- Ward functional was ex- 
tended also for superconductors,'^''^ Bose-Einstein condensates,'® and relativistic 
nonequihbrium phenomena in high-energy physics .'^1' •'33 

Let us introduce a nonequihbrium extension of the functional, i.e., ^ = ^[G*-'], 
in terms of S-matrix ()3-23p by 

^ = [ln(5c)]gkeleton,G'J(0)^G'J- = [(•^c)c - l] skeleton.G^iTO^G'i " ^^'^^^ 

Here, the subscript "skeleton" denotes retaining only skeleton diagrams without the 
self-energy insertion in the perturbation expansion for {•Sc)c — 1; and G^^^^> — )■ G*-? 
signifies replacing every G^^^^' by the renormalized one G^^ . See, e.g., Ref. lBip for the 
proof of the second equality that taking the logarithm of {Sc) amounts to retaining 
only connected diagrams of (Sc) — 1- Note ^* = ^. 

The Feynman rules for ^ are identical with those of Green's function above. It 
turns out that the exact self-energy is obtained from this <P by^ 

r«(l,2) = ±(-l).«;^^. (3.31) 

This functional differentiation graphically amounts to removing a Green's function 
line from every Feynman diagram for (p in all possible ways. The factor it originates 
from a reduction of a closed particle, and the factor (—1)*+-' reflects two fa's in Eq. 
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(|3^22]) . It follows from Eq. (IHaH]) . ^* = cp, and Eq. (f3^3B that Z" obeys 

^(l,2) = -firt(2,l)fi. (3-32) 

Now, our self-consistent perturbation expansion denotes: (i) carrying out the 
expansion of Eq. (|3-30p up to some order or retain some partial infinite series, and 
(ii) determining G and X! self-consistently using Dyson's equation (|3-20|) and relation 
()3-3ip . In this procedure, G is determined with E using Eq. ()3-20p . whereas IJ is 
given in terms of G as Eq. (j3-3ip : hence, the word "self-consistent". The advantages 
of the approach are summarized as follows. 

(a) It becomes exact when all the Feynman diagrams are retained, thus having the 
same structure as the exact theory. 

(b) It satisfies various conservation laws, as shown by Baym.'^ Indeed, it is clearly 
in the category of (^-derivable approximation. 

(c) Vertex corrections are taken into account automatically. 

(d) Two-particle and higher-order correlations can be calculated from any approxi- 
mate (P, i.e., providing some ^ amounts to determining the BBGKY hierarch}^^! 
completely. 

Point (b) makes the approximation quite suitable for describing dynamical phenom- 
ena. Indeed, no other systematic approximation scheme with the property seems to 
have been known to date within quantum field theory. It is worth pointing out that 
the Boltzmann equation can be derived by the approach with further approximations 
such as the quasiparticle approximation for the spectral function^ see ^5.11 below 
on this point. 

We now consider the first few series for ^ and write down the corresponding IJ 
explicitly. As already mentioned, Feynman rules for ^ are identical with those of 
Green's function given around Eq. ()3-26p . 

To begin with, we consider the first-order contribution to ^. Topologically dis- 
tinct diagrams are given in Fig. 5, where thick lines denote G; arrows, space-time 
arguments, and path indices have all been omitted. The corresponding analytic 
expression is given by 

2 

^(1) = ^1"^ I'di jdV {-lf-^V{l-l') [G**(l, 1)G^^(1', 1') ± G"(l, l')G"(l', 1)] . 

(3-33) 
It should be noted that there is an extra factor 1/2 in Eq. (j3-33p compared with Eq. 
(|3-27p . In general, the factor l/2n remains uncanceled in the nth-order contribution 
to ^ due to the fact that diagrams of <? are closed so that the 2n Green's functions 
are equivalent there'll 

The first-order self-energy is obtained from ^^^^ using Eq. ()3-3ip . Noting the 



Fig. 5. Topologically distinct first-order diagrams for <l>. 
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Fig. 6. Topologically distinct second-order diagrams for #. 



rule (d) below Eq. ()3-26p . we can put the result into the concise form: 

^«(l,2)=f3r«^(l,2), (3-34) 

with 

U^^{1, 2) = ±5(1, 2) / dl' V{l-l')ihG^\l', 1') + V{l-2)ihG^^{l, 2). (3-35) 



Except for the difference between G^"^^^' and G^^, Eq. (j3-34p coincides with Eq. 
(|3^ . as it should. 

We next consider the second-order contribution to (p. Topologically distinct 
diagrams are enumerated in Fig. 6. The corresponding analytic expression is given 
by 



^(2) 



(ih) 



2 2 



JJ=1 



^(-iy+-'' dl dl' d2 d2' V{1-1')V{2-2')[G'\1,2)G^'{2,1) 



xG'^ {!' ,2')G^\2' ,1') ± G'^ {1,2)G^\2,1')G'^ {!' ,2')G^\2' ,1)]. (3-36) 

The second-order self-energy is obtained from ^(^' using Eq. ()3-3ip as 

Z'*-'(2)(l,2) = {ihf I dl' f d2' V{l-l')Vi2-2')[±G'^ {1,2)G'^ {!' ,2')&\2' ,1') 

+G'^1, 2')G^\2', 1')G'^{1', 2)] . (3-37) 



Let us substitute Eqs. (|3-15cp and (|3-15dp into Eq. (|3-37p with i = j = 1 and compare 
the result with the expressions of i7^^'^^ and U'^^^'^' . We thereby obtain 

i;ii(2)(i,2) = (i/i)2 f dl' f d2' V{l-l')V{2-2'){e(h - t2)[±G^Hl,2)G^\l' ,2') 

xG^\2', 1') + G^\l, 2')G^\2', l')G^\l', 2)] + 0(t2 - h) [±G'^\l, 2) 
xG^^{l', 2')G^\2', 1') + G^^{1, 2')G^\2', l')G^^(l', 2)] 

= e{h - t2)u^^^^^ (1, 2) + eit2 - ti)i:i2(2) ^^^ 2). (3-38) 

Similarly, we have 

i:22(2)(i, 2) = e{t2 - ti)i72i(2)(i, 2) + e{ti - t2)u^^^^\i, 2). (3-39) 



Thus, we realize that the relations ()3-15cp and (|3T5dp between the diagonal and 
off-diagonal elements of G also hold among the elements of IJ^"^' . They are also 
satisfied generally among the elements of Z"'"' (n = 3, 4, • • • ), as may be checked by 
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inspecting higher-order contributions. Combining this fact with the first-order result 
(|3-34p . we conclude that the diagonal elements of IJ can be expressed generally as 

r"(l, 2) = r«^(l, 2) + e{t2-ti)S^\l, 2) + eiti-t2)U^\l, 2), (3-40a) 

^22(1, 2) = -U^^{1, 2) + e{ti-t2)U^\h 2) + e{t2-ti)S^H'^, 2). (3-40b) 



Hence, it follows that, besides Eq. (|3-32p . IJ also satisfies 

TrZ'(l,2) =TrfiZ'(l,2). 



(3-41) 



Thus, 17 carries the same symmetry as Eq. (|3-18|) of Green's function. Equation 
()3-40p tells us that we only need to calculate S^"^ and S'^^ of 17 for n > 2. 

3.5. Keldysh rotation 

To see the structures of nonequilibrium Green's function (j3T4p and Dyson's 
equation ()3-20p more clearly, it is sometimes convenient to perform the following 
"Keldysh rotation" to f^G^Bm 



G^ = LhGD 



qR qK 

G^ 



L 



V2 



1 -1 
1 1 



(3-42) 



By using Eqs. (j3-15p and ()3-17bp . the new elements G^' '^ can be expressed explicitly 
in terms of G^^ and G^^ as 



G^(l, 2) = e{t, - t2) [G^\l, 2) - G^\l, 2)] , 
G^(l, 2) = -e{t2 - ti)[G^\l, 2) - Gi2(l, 2)] , 



(3-43a) 
(3-43b) 



G^(l,2) = G21(1,2) + Gi2(l,2) = -'-{M'^)i,i{2)±i,],{2)^n{l)). (3-43c) 

The functions G and G thus introduced are the retarded and advanced Green's 
functions, respectively.- ' Using Eq. ()3-17ap . one may check that G ' ' satisfies 



[G^{1, 2)]* = G^(2, 1), [G^(l, 2)]* = -G^(2, 1 

which is expressed alternatively in terms of G^ as 

6^(1,2) =f2[G'^(2,l)]tf2. 
It is convenient to carry out the same transformation for E as 

I]^{l,2) = Lf:it{l,2)D -- 



17^(1,2) 17^(1,2) 
17^(1,2) 



Equation ()3-40p enables us to write down the elements of E as 

i7^(l, 2) = 17^^(1, 2) + e{ti - t2) [1721(1, 2) - E'\l, 2)] , 
17^(1, 2) = i7«^(l, 2) - e{t2 - ti) [1721(1, 2) - 17^2(1, 2)] , 
17^^(1,2) = 1712(1,2) + 1721(1,2). 



(3-44) 
(3-45) 

(3-46) 



(3-47a) 
(3-47b) 
(3-47c) 
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It follows from Eq. ()3-32p that IJ also satisfies 

rK(l,2) = f2[rK(2,l)]tf2. (3-48) 

With these preliminaries, we perform the Keldysh rotation to Dyson's equation 
(j3-2Up . Let us multiply Eq. (j3-2Up by Lfa and U from the left and right, respectively, 
and use LU = 1 with 1 the unit matrix. We thereby obtain Dyson's equation for G 
as 

ih-^ - kA G^(1, 2) - /" IJ^ih 3)G^(3, 2) d3 = 16{1, 2). (3-49) 

It follows from Eqs. (j3-42p and (|3-46p that: (i) the (2, 1) element of this equation 
vanishes identically, and (ii) the equation satisfies a symmetry relation corresponding 
to Eqs. (j3-45p and (j3-48p . Thus, the Keldysh rotation makes Dyson's equation more 
transparent. It should be noted at the same time that the perturbation expansion 
for G or IJ is performed more conveniently in terms of G''^ . 

The transformation (|3-42p adopted here is due to Larkin and Ovchinnikov^ and 
different in form from the original one by Keldysh.-' However, the physical contents 
are the same between the two transformations; moreover, the one by Larkin and 
Ovchinnikov has been used more frequently in condensed matter physics, especially 
in the field of superconductivity.'^ 

§4. Quantum transport equations and nonequilibrium entropy 

Dyson's equation (|3-20p and Eq. (j3-3ip for the self-energy form a closed set of 
self-consistent equations for Green's function G. Solving them, we can trace the 
nonequilibrium time evolution of the system. Over the last decade, active investiga- 
tions have been performed in the field of high-energy physics to solve those equations 
numerically for several approximate <^'s without any further approximations.'^T^ 

Here, we reduce those equations further to obtain quantum transport equations 
in the phase space based on the standard prescriptions of the Wigner transformation 
and a subsequent gradient expansion. We thereby remove some variables of Dyson's 
equations, which are irrelevant in many cases. The approximation was numerically 
checked to hold excellently over a wide range except for an initial time interval 
much shorter than the time scale for thermalization.'^ The resultant transport 
equation will be shown to contain information of entropy flow, i.e., it enables us to 
define nonequilibrium entropy that evolves with time so as to be compatible with 
an equilibrium expression. Thus, we will be able to clarify on what conditions the 
concept of entropy holds. The derivation below will be based on Ref. [25]) by Ivanov 
et al. andRef.lMD- 



4.1. Wigner representation 

The Wigner representation was introduced by Wigner in 1932 to consider quan- 
tum corrections to classical statistical mechanics .l^'l^ The representation enables 
us to extend the concept of "phase space' in classical statistical mechanics to quantum 
statistical mechanics. It also has a close connection with the Weyl quantizatioiP^ 
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and was used by Moyal in 1949 in his statistical formulation of quantum mechan- 
ics in phase spacePSJ Thus, the Wigner representation has had quite an impact in 
clarifying the foundation of quantum mechanics as well as its connection to classical 
mechanics. Moreover, it provides us with an indispensable tool for deriving quantum 
transport equations. 

To be specific, the Wigner transformation to G{1,2) is defined by 

G{ps,ri2ti2)= /'(i3fi2dti2G'(l,2)e-*(P-*^i2-rti2)/;i^ (4.1a) 

where ri2 = (ri + r2)/2, ti2 = (ti + t2)/2, fu = ri - r2, and ti2 = ti - t^. 
Equation ()4Tap is nothing but the Fourier transform of (5(1,2) in terms of the 
"relative" coordinates (ri2,ii2)- The quantity G{pe,ri2ti2) thus introduced is called 
the Wigner representation. The inverse transformation of Eq. ()4-lap is given by 

G(l, 2) = y ^^ G{pe, n2ti2) e^(P-i2--*"i2)/^ (4.1b) 

The function Wigner introduced in 1932 corresponds to the case with no time depen- 
dence.^^^ It is sometimes called the Wigner quasi-probability distribution with only 
(p, r) as arguments. Note that G{pe,rt) is different in form from (5(1,2); using the 
same symbol may not cause any confusion. Below every G given without arguments 
will denote G{pe,rt). 

Equation (j3-15p tells us that the independent components of G are G^^ and G^^ . 
It is useful for later purposes to express their Wigner representations in terms of two 
alternative functions. Let us first define the spectral function A by 

^(1,2) ^ i [G^\l,2) - G'\l,2)] = ^{i;nil)i'li2)T4'lm4'uil)), (4-2) 

which satisfies A*{1,2) = ^4(2,1), as shown with Eq. ()3T7ap . It follows from this 
symmetry and Eq. ()3-4ap at equal times that the Wigner representation of ^(1,2) 
has the following properties: 

f°° de 
A*{pe,rt) = A{pe,rt), / — A{pe,rt) = l. (4-3) 

J-oo 2vr 

Equation (|4-2p extends the equilibrium spectral function in the Matsubara formal- 
jgjjPUi ,13911 to nonequilibrium cases. We next introduce the distribution (j) directly in 
the Wigner representation as 

/-I 
d'n2-^(4(2)V^H(i)>e-'(p-^^^-^*^^)/'' 

npe,n2Ti2) - — A{pe,rutu) = <P* (pe , r^2t,2) . 

(4-4) 
As will be seen later, 4>{p£, rt) in equilibrium reduces to the Bose or Fermi distribu- 
tion function f{e) = {e^'^^l)~^. The two functions A and (j) form an independent 
pair of functions in G alternative to G^^ and G'^^. Thus, in nonequilibrium cases, we 
also need to determine the distribution function besides the spectral function. 
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In the conventional study of transport equations,—' one almost immediately 
integrates A out completely with some approximations such as the quasiparticle 
and quasiclassical approximations. One thereby obtains a simplified equation for (p 
alone. However, there may be cases where the approximations are not appropriate, 
such as those where the energy-level spacing is large or the density of states changes 
substantially. We proceed here by retaining A as it is, which also helps us to clarify 
the structure of the equations to be solved and its connection with the equilibrium 
formalism. 

The two quantities A and completely determine G. Indeed, it follows from 
Eqs. (f3T5]) . (|T2]1 . and (I33D that G^^ ^nd G^^ are expressed as 

G^^(pe, rt) = =FM(pe, rt)(j){pe, rt), (4-5a) 

G^^ (pe, rt) = -iA{pe, rt) [I ± 4>{pe, rt)] . (4-5b) 

Moreover, the Wigner representation of Eq. (|3-43p can also be written in terms of A 
and (p as 

G^{pe,rt) = [G^ipe,rt)]* = T ^ ^(P^^ (4.6,) 

G^{pe,rt) = -iA{pe,rt)[l ± 20(pe,rt)], (4-6b) 

with e+ = e + iO-|-. Note that Eq. ()4-6ap is a convolution of A and e^ , where e^ is 
the Fourier transform of step function (j3-16p . Thus, G^ contains information on the 
distribution function, whereas G^' can be expressed solely in terms of the spectral 
function. 

The same transformation is also possible for the self-energy matrix IJ, whose 
diagonal elements can be expressed in terms of the off-diagonal elements as Eq. 
([Mo]) . Let us write T^^ ^^^^ ^21 ^^ 

U'^'^ipe, rt) = ^fiA^ipe, rt)(j)s{pe, rt), (4-7a) 

E'^\p£,rt) = -iAjj{p£,rt)[l ± (f)Eip£,rt)], (4-7b) 

where A^ and (pz; are two alternative independent functions. It then follows from 
Eqs. (|3-35p and (|3-47p that the elements of Z"^ are expressed as 

E^ipe, rt) = [uHpe, rt)]* = E^^ip, rt) + T^ ^P^^ (4.8a) 

E^{p£,rt) = -iAEip£,rt)[l ±2(l)z{pe,rt)]. (4-8b) 

Finally, it is worth pointing out that the operator: 

G^\l,3) ^ (^ih^^- k,y (1,3), (4-9) 

where Ki is defined by Eq. ()3-3p . can also be transformed into the Wigner represen- 
tation as 

G^\pe,rt)^e-^-U{rt)+f,. (4-10) 

Equation (|4-10p is useful when converting Eq. (|3-49p into the Wigner representation. 
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4.2. Groenewold-Moyal product and gradient expansion 

We have introduced the Wigner representation as Eq. ()4-lap . where G{1,2) can 
be regarded as a matrix with the arguments 1 and 2. With this viewpoint, the 
self-energy term in Dyson's equation (|3-49p is a matrix product. Hence, we need 
to clarify how a matrix product changes through the Wigner transformation before 
expressing Eq. (j3-49p in the Wigner representation. The answer to this question is 
given in terms of two arbitrary matrices C(l, 2) and D{1, 2) by 

I C{1, 3)Z)(3, 2)d3 = I ^^C{pe, mtu) * D{pe, rutu) e^(P-l2--*l^)/^ (4-11) 

Here, the operator * is defined in terms of dr = d/dr and dt = d/dt by 

C{pe,rt) * D{pe,rt) = C{pe,rt) exp — ( drOp — d tO e — dpOr + d sOi 

xD{pe,rt), (4-12) 

where the left (right) arrow on each differential operator denotes that it works on 
the left (right) function. Equation ()4-lip was derived independently by Groenewold 
in 194#9 and by Moyal in 1949.1221 We call it the Groenewold-Moyal product. 

By considering only the space coordinates for simplicity, Eq. (|4-lip is proved 
as follows. Let us express C(ri,r3) and D{rs,r2) on the left-hand side as Eq. 
(|4-lbp . We next expand the "center-of-mass" coordinates ri3 = T12 + 2^32 and r32 = 
ri2 - ^ns in Taylor series from ru as C{p,ri2 + ^f32) = exp[ir32 • dri2]C{P,ri2) 
and D{p',ri2 — ^rn) = exp[— iri3 • dr 12] D{p',r 12). We then remove the relative 
coordinates f32 and ri3 in the expansions as e^^''^^^'^^ exp^—^fi^ ■ S^^j] L>(p', ri2) = 
exp [y9p- 9^.12] e*^'^^^'^-D(p',ri2), and subsequently perform partial integrations over 
p and p' as C7(p, ri2) exp [f Sp-a^.J e*P*^"A ^ C7(p^ ^^2) exp[-f ^p-a^.J e^P-*^"//i, 

Finally, we carry out the integration over r^ with only plane waves in the integrand 
to obtain the factor 6pp' e^P''^^'^ . 

Retaining up to the first order in the gradient expansion of Eq. (|4-12p yields 

C*D^CD+'-^{C,D}, (4-13) 

where {C, D} denotes the generalized Poisson bracket: 

srn\ = — —- — — -— — ^^ 

dr dp dt de dp dr de dt 

This approximation will be called the first-order gradient expansion. It will hold 
excellently when the microscopic length scale /m and time scale im of the system 
are much smaller than the macroscopic length /m and time tyi characterizing the 
inhomogeneity of the system. For a system of low-temperature fermions with the 
Fermi wave length k-p and Fermi energy ep, for example, the microscopic length 
and time are given by Im ~ l/^F and tm ~ ^/epj and those in a classical gas 
are the mean spacing between two particles and h/k^T, respectively. When the 
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space-time variations of the system occur over scales much longer than them, Eq. 
(j4-13p will hold excellently. This condition is certainly met over a wide range of 
nonequilibrium systems. Indeed, the approximation was numerically checked to hold 
over a wide range except for an initial time interval much shorter than the time scale 
for thermalization.--' 

4.3. Quantum transport equations 

With these preliminaries, we now derive quantum transport equations from 
Dyson's equation ()3-49p . It is useful for this purpose to write its first term as a 
matrix product with Eq. (|4-9|) . Then, we can use Eq. (j4-ll|) to obtain the Wigner 
representation of Eq. (|3-49p as 

{G^^i-IJ^)*G^ = 1, (4-15) 

where Gq^ is defined by Eq. (I¥T0|) . 

We next adopt approximation ()4-13p for Eq. (|4T5p . By noting Eqs. ()3-42p and 
()3-46p . the (1,1) element of the resultant equation can be written down explicitly as 

(Gq-i - i:^)GR + y{Go-i - r^, G^} = 1. (4-16a) 

The (2,2) element is obtained from the above with the replacement R— t-A in the 
superscripts. Taking its complex conjugate subsequently with G = G^ in mind, 
we have 

(G^i - i;^)GR - y{Go-^ - r^, G^} = 1. (4-16b) 

This equation can also be derived directly from the (1,1) element of the right-hand 
Dyson's equation: G^ * {Gq 1 — IJ^) = 1. Let us add the above two equations. We 
then obtain 

G^ = {Gq^ - U^)-\ (4-17) 

Noting Eq. ()4-6ap . we realize that this equation determines the spectral function A 
for a given S^. 

Equation (|4-17p has been derived by adding Eqs. (|4-16ap and ()4-16bp . However, 
the (1,1) and (2,2) elements of original Eq. ()3-49p are complex conjugates to each 
other so that they are equivalent mathematically. Hence, one may ask whether this 
equivalence is still kept between Eqs. ()4-16ap and (|4T6bp obtained with the first- 
order gradient expansion. To answer this question raised by Ivanov et al.,'^ we 
alternatively subtract Eq. ()4-16bp from Eq. (|4-16ap and subsequently substitute Eq. 
(I¥T7)1 . It yields 

= {(G^)-\G^} = -(G^)-2{G^,G^}, 

which is just the identity {G^,G^} = 0. Thus, Eqs. (|4-16ap and ()4-16bp have 
been checked to be identical also within the first-order gradient expansion. This 
consideration shows that G^ from Eq. (|4T7p is correct up to the first order. 

We next focus on the (1,2) element of Eq. (l¥T5]l . Noting Eqs. (l3^i2D and (f3^i6D 
and adopting first-order gradient expansion ()4-13p . we obtain 

(G^i - i:^)GK - r^G^ + y {Gq-i - E^, G^} - ^{i:^, G^} = 0. (4-18a) 



Introduction to Nonequilibrium. Statistical Mechanics 25 

Its complex conjugate can be expressed with Eqs. ()4-6p and ()4-8p as 

- (G^i - U^)G'^ + E'^G'^ + -{G^i - E^, G^} - -{s"^^ G^} = 0. (4-18b) 

Equation ()4-18bp may also be derived directly from the (1,2) element of the right- 
hand Dyson's equation G^ * {G^^l - S^) = 1. 

Let us add Eqs. (|4-18ap and (|4-18bp and make use of Im{G(^^ - 17^, G^} = as 
well as Eqs. (|4-6p and (|4-8p . We thereby obtain 



;o'i - Rer^, ^0} - {^i:<^i:, ReG^} = ^IM^ 



The left-hand side of this equation consists only of space-time derivatives, which are 
first-order in the gradient expansion. Hence, it follows that the term (/)^ — (p on 
the right-hand side is also of first-order; as seen below in ^4.61 the right-hand side 
may be identified as the collision integral in the transport theory, which vanishes in 
equilibrium. Hence, it follows that we may replace (px: by (p on the left-hand side as 
adopted by Botermans and Malfliet.^^''®^' The approximation yields 

{Gq^ -Rer^, Acj)} - {^170, ReG^} = C, (4-19) 

where C is the collision integral defined by 

Equation ()4T9p determines the distribution function <p for given A and S. 

Alternatively, one may subtract Eq. (|4-18bp from Eq. ()4-18ap . Then, a procedure 
similar to that for deriving Eq. (j4T9p yields 

{Ae(P,A} - {Ae,AcP} _ AscPe^^G^ - {G^^ -KeE'^)A<p 

4 " h ■ ^^'^^^ 

This equation is nothing but Eq. ()4-19p in disguise. Indeed, multiplying Eq. ()4-2ip 
by Ae / {Gq'^ -KeS^) reproduces Eq. (j¥T9]l . This can be checked with Eq. (f¥T71l by 
expressing ReG^ and A = -2ImG^ in terms of M^G^^-ReZ"^ and ^i: = -2ImZ'^ 
as 

RprR = ^^ A = -^^ 

M2 + (^^/2)2' ^ M2 + (^^/2)2' 

and substituting them into the two apparently different equations. Thus, we have 
seen that the equivalence between Eqs. ()4-18ap and (|4T8bp is recovered appropriately 
with the replacement (ps^-fj) i^i the space-time derivatives. 

Equations ()4-17p and ()4T9p form a closed set of equations for the two unknown 
functions A and (p. Moreover, the replacement (p — )• /(e) = (^e'''^ =F 1)""*^ in those 
equations transforms: (i) Eq. (|4-17p into Dyson's equation for the retarded Green's 
function in equilibrium; and (ii) Eq. ()4-19p into the trivial relation = 0, i.e., both the 
space-time derivatives and collision integral become null with the replacement (see 
^4.6l below). Thus, Eqs. (|4-17p and (|4-19p form a natural extension of the equilibrium 
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Dyson's equation into nonequilibrium cases, where we also have to determine the 
distribution function (j) with Eq. (|4-19p . 

Finally, it is worth pointing out that a Wigner representation may not be positive 
definite in general. However, the following argument suggests that the spectral 
function A and distribution function <f> obtained using Eqs. (j4-17p and (j4-19p are 
both positive. First, they are real, as Eqs. ()4-3p and ()4-4p show. Second, it follows 
from the retarded nature of G (1,2) in Eq. ()3-43ap that all the singularities of 
G^{pe,rt) in Eq. (j4-17p lie on the lower half of the complex e plane. This implies 
that Im.U^{pe, rt) < and, hence, ImG^(pe, rt) < 0. By using Eq. ()4-6ap . the latter 
condition can be written alternatively as 

A{pe,rt)>0. (4-22) 

As for (p, it reduces in equilibrium to the Bose/Fermi distribution function, which 
is definitely positive. Noting that it obeys Eq. ()4-19p of the first-order gradient 
expansion so that its deviation from the equilibrium form will not be substantial, (p 
is also expected to be positive. 

4.4. Approximation for the self-energies 

We now discuss how to calculate the self-energies within the first-order gradient 
expansion of the self-consistent (^-derivable approximation. In brief, they should 
be estimated with the local approximation of neglecting the space-time derivatives. 
This may be realized by noting: (i) Eq. ()4-17p is correct up to the first order in the 
gradient expansion even with S^ of the local approximation, as pointed out in the 
paragraph just below Eq. (|4T7p : (ii) Eq. (|4T9p . which is composed of space-time 
derivatives and the collision integral, is first order by itself so that we should make 
use of the zeroth-order self-energies in the equation. 

We now write down the self-energies in the local approximation up to the second 
order in the perturbation expansion. Using Eq. (|3-5p . we can transform Eq. ()3-35p 
in the Hartree-Fock approximation into the Wigner representation as 

i7HF(p,rt) = ±ihj |^(^o ± Vp^p,)G^Hp'e',rt). (4-23) 

Similarly, second-order self-energy (|3-37p may be expressed with Eq. (|4-lip in the 
Wigner representation. Adopting the local approximation subsequently, we obtain 



U'^^''\pe,rt) = ^{hfJ 



{2nhY 



-\Vp-ps, ± Vp-p^\^{2-Khf5{p+p2-p-i-PA) 

k=2 " ' J 

x5{e+e2-e^-eA)G^\p2e2,rt)G'^{p^e^,rt)G'\pAeA,rt). (4-24) 

By writing G^'^ = ^iG<, G^^ = -iG>, U'^'^ = ^iU<, and U'^'^ = -iU>, this expression 
acquires exactly the same expression as Eq. (4-16) of Kadanoff and Baym.'^ 

Having given the self-energies as a functional of G^^ in the Wigner representation, 
we can solve Eqs. (j4-17p and (|4T9p self-consistently to trace their time evolutions. 
The exact procedure is summarized as follows: (i) With <^(pe, rt) given at a certain 
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time t, we calculate the spectral function A{p£, rt) self-consistently using Eq. ()4-17p 
and S^ = U^{p£,rt]A,4>) of the local ^-derivable approximation, (ii) The distri- 
bution function (j) at t + dt is determined subsequently with Eq. (|4-19|) . Hence, it 
follows that the explicit time dependence in the procedure originates from (j). 

4.5. Expression of entropy density 

Equation ()4T9p . which describes the time evolution of the distribution function, 
also contains information on the heat flow. Indeed, it enables us to derive an ex- 
pression of nonequilibrium entropy density and its continuity equation. To see this, 
let us multiply Eq. (|4-19|) by /ifeB(27r/i)~^ln[(l it (/>)/(/)], carry out an integration over 
pe, and make use ofIm.{GQ^ — I!^,G^} = and ln[(l ± (l))/(f)]d(p = da with 

a[(l)] = -(p\n(p±{l±(t>)ln{l±(t>). (4-25) 

We thereby obtain 

ds . dscou 



where s = s{rt), js=js{ft), and dscoii/dt are defined by 



d^pde 
d^pde 



^ ° o + Ae^ 

oe oe 

d{GQ^-ReS^) dReG^ 

dp dp 



(4-27) 

(4-28) 



^'^""-^^b/I^CIuI^. (4.29) 



dt J {2-Kh) 

Equations ()4-27p and ()4-28p may be identified as entropy density and entropy flux 
density. Indeed, if we put 4> — )• f{e) = (e"^'*"'^^ =F 1)~^ in Eq. (j4-27p and perform an 
integration over r, we reproduce an expression of equilibrium entropy derived in Ref. 
^2]) : see also Appendix [D] for the equilibrium expression. On the other hand, Eq. 
(|4-29p denotes entropy production per unit time and unit volume due to collisions; see 
the next subsection on this point. Thus, we have been able to define nonequilibrium 
entropy density so as to be compatible with equilibrium statistical mechanics. As 
will be shown in ^5.H it also embraces Boltzmann's nonequilibrium entropy density 
as the dilute high-temperature limit. 

Expression ()4-27p is different from the one obtained earlier by Ivanov et alW^ 
which contains extra terms stemming from space-time derivatives in the self-energies, 
i.e., terms due to memory effects. The corresponding expression of nonequilib- 
rium entropy density does not reduce adequately to the entropy density in equi- 
librium^^ ''^^' as (/» — )• /(e) = (e"^' bT^ zp 1)"-^. It was derived so as to be compatible 
with the equilibrium entropy by Carneiro and Pethick.'^ However, the derivation 
of Carneiro and Pethick was based on the zero-temperature perturbation expansion 
technique of Goldstone^^^ and suffers from an inappropriate treatment of energy de- 
nominators, as pointed out in Ref. I^6|) . Thus, it is Eq. (j4-27p that is compatible with 
equilibrium statistical mechanics. 
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4.6. H-theorem and thermodynamic equilibrium, 

Here, we show that Eq. ()4-26p satisfies the i?-theorem, i.e., the law of increase 
in entropy, within the second-order perturbation expansion. 

Let us substitute Eq. (|4-24|) into Eq. (|4-2Up and the corresponding expression 
of C subsequently into Eq. ()4-29p . Using Eq. ()4-5p , we further express the result in 
terms of A and (j) and symmetrize the formula with respect to integration variables. 
The second-order quantity ds).jyjdt is thereby transformed into 



dt 



4 r 

X^lA2^3^4[(l±</'l)(l ± 02)</'3</'4-01</'2(l ± h){l ± </'4)] 
"^%,^2(1±^3)(1±04)' ^'-'^^ 

with Aj = A{pjej,rt), etc. Noting (x — y) ln{x/y) > 0, which holds for any positive 
X and y, we conclude ds^J^^/dt > 0. This is nothing but a quantum- mechanical 
extension of Boltzmann's i^-theorem.'^^''^ The equality in Eq. (|4-3Up holds when 

01 02 03 04 

is satisfied. Noting the delta functions in Eq. (|4-30p . we realize that ln[(lib0i)/0i] 
must have linear pi and ei dependences as 

lni4^ = a + /3(ei--j;-pi), (4-31) 

01 

where a, /3 and v are arbitrary functions of rt. We thereby obtain, for = 0^'*^) that 
makes the collision integral vanish, the expression: 

0('^)(pe,rt) = -r— r^ , ^, .^ , ^, , (4-32) 

which is exactly the local equilibrium distribution. 

Thus, we have proved the //-theorem within the second-order perturbation and 
also provided an expression for the distribution function 0^'*^) that makes entropy 
production vanish. It can also be shown easily within the second-order perturbation 
that this 0*-^''^ makes the collision integral (j4-20p vanish. 

It still remains to be clarified whether dscoii/dt > does hold up to the infinite 
order. However, it seems quite natural physically to expect it. 

§5. Quasiparticle and quasiclassical approximations 

In many cases, quantum transport equations (|4-17p and (|4-19p are further sim- 
plified by integrating out the spectral function A completely. There are a couple 
of standard approximations to carry it out, i.e., the quasiparticle and quasiclassi- 
cal approximations. We explain them here in detail using second-order self-energy 
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()4-24p . It is thereby shown that quantum transport equations embrace the Boltz- 
mann equation as a Hmit. It will also become clear under what conditions the concept 
of "distribution function in the phase space", which is frequently assumed to exist 
from the very beginning in the classical transport theory}^ is justified. 

5.1. Quasiparticle approximation 

This approximation becomes excellent whenever there is a sharp (5-function- 
like peak in the spectral function A. We first consider the weak-coupling limit as 
a typical example. Here, we may drop the self-energy in Eq. (|4-17p . Hence, the 
spectral function A = — 2ImG^ is obtained from Eq. (|4-1U|) as 

A{pe, rt) « 2tt6 {e - Cp) , ip{rt) = ^ + U{rt) - /x. (5-1) 

As for Eq. ()4-19p . we can neglect terms with S on its left-hand side to the same 
order of approximation. By noting Eqs. (|4-10p and (j4-14p . the resultant equation 
can be written down explicitly as 

d{A<^) ^ p d{Acp) ^ dU d{A(t>) dU d{A(t>) ^^ 
dt m dr dt de dr dp 

We next divide Eq. (j5-2p by 27r and perform an integration over — oo < e < oo. 
The corresponding equation can be expressed concisely in terms of the distribution 
function in the phase space: 

/■°° de 

f{p,r,t)= —A{p£,rt)(t){p£,rt) = (t){pip,rt), (5-3) 

J-oo 2vr 

^ + ^ • ^ - — • ^ = / [/I (5-4a) 

dt m dr dr dp ^ 

The collision integral Ip[f] can be written down by substituting Eqs. (j4-5p . ()4-24p . 
and (|5-ip into Eq. (j4-2Up and carrying out the integration over e as 

iv\f] = Y n y (^i^p-P3 ± vp.pf{2^)H{p +P2-P3- P4) 

xSiCp + ip, - ip, - ip,) [(1 ± /)(! ± f2)hh - //2(1 ± /3)(1 ± U)] , (5-4b) 

with fj = f{pj,r,t). Equation (|5-4p may be regarded as the Boltzmann equa- 
tion with quantum effects. Indeed, if we neglect quantum effects in Eq. (|5-4bp as 
±Vp-p^V*_p^ — )• and (1 ± /) — )• 1, we obtain the standard Boltzmann equation 

in the dilute classical limit .1^ Note that the function / here is dimensionless and 
normalized so as to reproduce {e^^^ =F 1)~^ in equilibrium. 

The corresponding expressions of entropy density, entropy flux density, and en- 
tropy production are obtained by adopting the same approximation in Eqs. (|4-27p . 
(j4-28p , and ()4-30p , respectively. They are given in terms of Eq. (|4-25p by 

/dp 



30 



js = k 
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{27rh) 



3 -^[fh 



(5-6) 



ds 



(2) 
coll 



ksh^ 



dt 



n 



fj T) ' 

—^\Vp,-p,±Vp,^pf{27r)^6{pi+p2-p3-p4) 



xSi^p + ^P2 - ^Pa - ep4) [(1 ± /l)(l ± /2)/3/4 - /l/2(l ± /3)(1 ± /4)] 
^j^(l±/l)(l±/2)/3/4 



/l/2(l±/3)(l±/4)' 



(5-7) 



Equation ()5-5p is the familiar entropy density for noninteracting systems!^ These 
expressions can be obtained alternatively by multiplying Eq. ()5-4p by k-Q ln[(l it /)//] 
and performing an integration over p. Boltzmann's entropy and //-theorem are 
reproduced from the above with ±.Vp-p^Vp_p^ — )• and 1 it / — )• 1. 

Thus, we have seen that a transport equation for the distribution function / 
in the phase space results naturally as the weak-coupling limit of quantum trans- 
port equations (|4-17p and (|4-19p . It should also be noted that the weakness of the 
interaction is not essential for the quasiparticle approximation. To be specific, it 
is sufficient for the approximation to hold that the imaginary part of S^, i.e., the 
lifetime of quasiparticles responsible for the excitation, is sufficiently small in the 
relevant energy range. Thus, the quasiparticle approximation can also be applied 
excellently for normal Fermi liquids at low temperatures.'^lJ' We will discuss it in 
detail below. 

First, let us neglect terms with Ax; = — 2ImZ'^ on the left-hand side of Eq. 
(|4-19p . We also expand Gq — ReZ'^ up to the first order from e = and p = Pf 
{pp: Fermi momentum) as 



Gq^ - ReU^ 



e-?P 



Here, a and ^p are the renormalization constant and quasiparticle energy: 
1 . dReU^ 



1 



ip = a 



de 



£ = 0,P=Pf 



^ + [/ + Rei7^(0, 0) - /i 
2m 



+ VY- (p-Pf), 



respectively, with the Fermi velocity vp defined by 



IPF dReU^ 

VF = al \ 

\ m op 



(5-8) 

(5-9a) 
(5-9b) 

(5-9c) 



£ = 0,P=Pf / 



The spectral function A = — 2ImG in this approximation is obtained from Eq. 
(|4-17p as A{pe,rt) ss 2TTa5 {e — ^p). Let us introduce the distribution function / as 



f°° de 
f{p,r,t)= - — A{pe,rt)(j){pe,rt). 

J-oo 2vra 



(5-10) 
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We then obtain Eq. ()5-4p from Eq. ()4-19p using exactly the same procedure as above 
for the weak interaction. However, the factor Vp-p^ it Vp-p^ in Eq. ()5-4bp is to be 
replaced by the effective interaction between quasiparticles, as described in detail in 
Ref. [3T]) . for example. 

5.2. Quasiclassical approximation 

For Fermi systems at low temperatures, there is an alternative approximation 
called "quasiclassical approximation" . It holds excellently even when the lifetime of 
quasiparticles is substantial, as for electrons in metals with strong impurity poten- 
tials or electron-phonon interactions. The approximation was introduced by Prange 
and Kadanoff for the case of electron-phonon interactions,'^ and extended sub- 
sequently to describe superconductors and superfluids in equilibriuni^''® and in 
nonequilibrium .'2SI''E5|i-168li jj^ ^j^g ^.g^gg q£ conventional superconductors, for example, 

the impurity scattering plays a crucial role in changing their properties in magnetic 
fields from type I to type 11,'^ and the lifetime may not be negligible in the lat- 
ter cases. It was Eilenberger*'^- who first adopted the approximation to integrate 
out an irrelevant variable from the Gor'kov equations in equilibrium. He thereby 
obtained quasiclassical equations of superconductivity, i.e., Eilenberger equations, 
which extend the Ginzburg-Landau equations™-* over whole temperatures and mag- 
netic fields. Those equations have been used extensively to clarify and understand 
the properties of inhomogeneous superconductors and superfluids.'^''^ 

In the quasiclassical approximation, one carries out an integration over ^p defined 
by Eq. ()5-9bp instead of an integration over e in the quasiparticle approximation. 
To be specific, proceed here as follows. Let us write down Eq. ()4T9p explicitly with 
Eqs. (|TOll . d^BD, and (l5^ as 

ld{A(t)) v^ d{A^) d{U+RfiU^)d{A(^) d{U+ReS^) d{A(t)) 
a dt a dr dt de dr dp 

+i{r^2 ReG^} = C, (5-11) 

where we have approximated A^cf) ~ —iH^"^ on the left-hand side. The term 
dIJ{p£, rt)/dp is expected to be of the order of \IJ{pe, rt)\/pp so that we may neglect 
the IpI dependence in IJ as 

E{pe, rt) f« IJippe, rt). (5-12) 

By noting luiU = —2 Ax; and Eq. ()5-8p . the corresponding spectral function A = 
— 2ImG^ is obtained as 

a'^Ax;{pF£,rt 

{e-^p)^+ [^aAs{pFe,rt)] 

The IpI dependence in Eq. ()5-13p lies only in ^p and may be approximated excel- 
lently using a Lorentzian in terms of ^p with the width i^'^Aj] and the area 27ra. 
With this observation, we introduce the distribution function in the quasiclassical 
approximation as 

P''{pFe,rt)= lim / --^A{p£,rt)(t){pe,rt). (5-14) 

ic^ooj_^^ 27ra 



Mpe, rt) « , — ^: '::'^'r:' ,..2 - (5-13) 
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Thus, the distribution function here has e as an argument instead of |p| in the 
quasiparticle approximation, and the momentum has only the angular dependence 
on the Fermi surface. We next perform an integration over ^p in Eq. (j5-lip in the 
same way as in Eq. (|5-14p . It then follows that the fourth term on the left-hand 
side vanishes upon the integration. Moreover, the ^p integration of ReG^ may be 
estimated with Eqs. ()4-6ap and ()5-13p as 

lim /"^' ^ReG^(pe,rt) =P H —^—. (5-15) 

Cc^oo J_^^ 27ra 7_oo 27r e - e' 

Hence, it follows that the fifth term with derivatives of ReG^ vanishes after the ^p 
integration. As for the collision integral C obtained by substituting Eqs. ()4-5p and 
(l¥^ into Eq. (I¥^ . let us introduce the density of states by N{(,) = f j£^S{C - 
^p) and subsequently approximate it as N{^) ~ N(0) to carry out the ^ integration. 
We thereby obtain the quasiclassical transport equation as 



with 



I^V^l - y[iV(0)]3(2vra)^n [ ^^K.-P^.-^^ 



i=2- 

X5{PF +PF2 -PF3 -pF4)'5(e + e2 -£3 -£4) 

X [(1- mi-fr)frfr-r'fri^-fr)i^-fr)], (5-17) 

where df2j denotes an infinitesimal solid angle along ppj . 

Thus, we have obtained a transport equation in the quasiclassical approximation. 
As already pointed out at the beginning of this subsection, it has an advantage of 
being applicable even for cases with short lifetimes. The derivation here has been 
carried out with second-order self-energy (|4-24p of a weak interaction. However, the 
approximation is justified even for low temperature Fermi liquids with considerable 
interactions, and also applicable to superconductors and superfluids. See Ref. [66|) 
for details. 

§6. Transport equations for electrons in electromagnetic fields 

We next focus on electrons in electromagnetic fields, where special consideration 
of the gauge invariance is necessary for deriving quantum transport equations appro- 
priately. To be specific, we should modify Eq. (j4Tp to introduce a gauge-covariant 
Wigner transformation in terms of the "center-of-mass" coordinate and make use 
of the corresponding Groenewold-Moyal product. We will discuss it below following 
Ref.ESj). 

We consider the following one-particle operator in place of Eq. ()3-3p : 



2m 



e 



-ihVi--A{l)\ +eAAl)-^i, (6-1) 

c J 
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where A(l) and ^4(1) are the vector and scalar potentials of electroniagnetism, 
respectively. The corresponding Dyson's equation is still given by Eq. (|3-49p : it is 
invariant through the gauge transformation: 

A(l) ^ A(l) + Vix(l), ^4(1) ^ ^4(1) - -^, (6-2a) 

c oti 

Gil, 2) ^ exp{i-^[x(l) - xm} G(l, 2), (6-2b) 

with e < and x(l) denoting an arbitrary scalar function. 

To discuss transport phenomena of electrons in solids under electromagnetic 
fields, we also need to incorporate interactions with phonons, impurities, etc., to 
obtain a finite conductivity. However, Dyson's equation for the electronic part is 
still given by Eq. (j3-49p with modified self-energies. We will proceed by assuming 
that such effects have already been taken into account into the self-energy matrix of 
Eq. ()3-49p . See also Ref . [691) on how to incorporate the periodic potential. 

Following the standard prescription to derive quantum transport equations, we 
first transform Eq. ()3-49p into a Wigner representation. If we adopt Eq. ()4Tp for this 
purpose, however, the resultant equation suffers from the flaws that: (i) the gauge 
invariance of the equation is not retained adequately; (ii) the magnetic Lorentz force 
of deflecting electrons is absent from the equation. This is because the transformation 
(|4-ip breaks the gauge covariance with respect to the "center-of-mass" coordinate. 

To remove these drawbacks, we adopt the "gauge- invariant Wigner transforma- 
tion" introduced by Stratonovich in 1956.^^' It is defined in place of Eq. ()4-ip by 



G(l, 2) = e'''^'''> J -j^ G{pe, r^tu) e^P-^---^'-^/'', (6-3b) 

where 1(1,2) is given in terms of the four vectors A = (A, —cA/^) and fi = (ri,ti) 
by 

I{l,2) = ^^P A{s)-ds, (6-4) 

with s denoting the straight-line path from r2 to fi . Equation ()6-3p was also derived 
by Fujita elegantly on the basis of the Weyl transformation.'^ 

The role of the extra phase factor in Eq. (]6-2bp may be realized from its change 
under the gauge transformation (|6-2ap : 

ei/(i,2) _^ exp{i^[x(l) - x(2)]} e''^'^'\ (6-5) 

Thus, it transforms in exactly the same way as G in Eq. (|6-2bp . thereby making 
(5(1, 2)e~'^(^'^' in Eq. (j6-3ap gauge invariant. It is worth looking at Eq. (j6-3p from 
a more general viewpoint for its extension to superconductors. Setting 1 = 2 in 
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3 




1 2 

Fig. 7. Closed contour C123. 

Eq. ()6-2bp tells us that (7(1,2) is gauge invariant with respect to the "center-of- 
mass" coordinate. Thus, the gauge dependence of (5(1,2) lies solely in the relative 
coordinate, which should be removed before the Wigner transformation. The factor 
g-«^(i,2) Qf (5(1, 2)e~*-'^(^'^) exactly performs this task. It has been showri^S that 
this idea can be applied to superconductors so as to describe appropriately that the 
pair potential behaves like an effective wave function of charge 2e with the gauge 
covariance in terms of the center-of-mass coordinate; see the last paragraph of this 
section for more details. 

We now modify the Groenewold-Moyal product ()4Tip so as to be compatible 
with Eq. (|6-3p . We carry it out here for the simplest case of static electromagnetic 
fields; see Ref. [68]) for a more general treatment. We first express (7(1, 3) and -D(3, 2) 
on the left-hand side of Eq. (j4Tip as Eq. ()6-3bp . We then have the extra phase factor 
g«/(i,3)+i/(3,2) jj-^ ^YiQ integrand on its right-hand side. Let us express it as a product 
of e^^^'^' and the other factor as 

g^/(l,3)+./(3,2) ^ ^i^l,2)^^4,^2^^ ^^^^ ^^I X{s) ■ ds, (6-6) 

where the contour C123 is given in Fig. 7. Using Stokes' theorem, we can express (/)i23 
in terms of the electric field E = — VA4 — -dtA and the magnetic field S = V x A 

as 

0123 = wr- [-cE ■ (ri3t32 - tisr32) - B ■ (ri3 x 7=32)] . (6-7) 

The extra phase factor in Eq. ()6-6p modifies Eq. (j4-lip into 

j C7(l, 3)Z?(3, 2) ci3 = e^^(i'2) j ^^ C(pe, ^2*12) * D{pe, rutu) e^(P-l2--*l^)/^ 

(6-8) 
where the Groenewold-Moyal product is now defined by 

C{p£,rt) * D{pe,rt) = C{p£,rt)exp — (dr-Op- dtOe- dpdr+ d eO t 

xexpjy -e^-(^p^,-K'^p)+^S-(Kx^p) I 
xD{pe,rt). (6-9) 



This modified Groenewold-Moyal product was derived independently in Refs. [68j) 
and 1 73 p . Its expansion up to the first order in the differential operators can also be 
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expressed as Eq. ()4-13p with the generahzed Poisson bracket: 

srn\ = — ^_^^_^ ^ dCdP 

dr dp dt de dp dr de dt 

(d£dD_dCdD\ e f dC_ dD\ 
\dp de de dp J c \dp dp J 

With these prehminaries, we now express Dyson's equation (j3-49p with Eq. (j6-l|) 
in the Wigner representation. Let us substitute Wigner representation (|6-3bp into 
Eq. (j3-49p . We then have space-time derivatives of 7(1,2). To estimate them, we 
express the straight-Hne path s in Eq. (|6-4p in terms of R= 2(^ + ^2) and r = fi — r2 
as s = R + rjf {—-^ < ?? < 2)- ^^ ^^^° expand A{s) in a Taylor series from s = R 
and subsequently perform the integration over rj order by order. We thereby obtain 
an alternative expression of 1(1, 2) as 

1(1, 2) = ^5(r • d^) f- A{R), g{x) ^ "-^^^^. (6-11) 

We next operate df^ = dff+ ^dj^ on 1(1,2) and make use of the identities r^d^Ay = 
ryduA^ + ry{d^Ay-duA^) (^u, 1/ = 1, 2, 3,4) mid g{x) + [g'{x) + \g{x)]x = e"'^ . The 
derivative 9/(l,2)/9ri^ is thereby transformed as 



a/(l,2) _ e 



dr 



l^i 



he 



^m(I) 



g'{r-d^) + -g{r-d^) 



E^ 



dA,{R) dA^{R) 
dRa dR^ 



A{l) + -{cti2E + fi2xB) 



Ai{l) + -ri2 • E 



:/x = 1,2,3 
:// = 4 



(6-12) 



Here, the first equality is an identity, whereas the second one holds only for static 
electromagnetic fields. 

We now substitute Eqs. ()6-ip and (|6-3bp into the first term in Eq. ()3-49p and 
follow the procedure of deriving Eq. ()4-lip . Comparing the resultant expressions 
with the definition of the * product in Eq. (|6-9p . we obtain the simple expressions: 

^^ p * G^(p., n2ti2) e^(p-^— *"-)/^ 



-ihVi--Ail) G^(l,2) 
c J 



=i^(l,2) 



d 
ih-- - eAi{l) 

Ot\ 



G^(l,2) = e 



i/{l,2) 



{2-KhY 

d^pde 



e*G^{pe,ri2ti2)e 



(6-13a) 

i(p-fi2~eti2)/fi 



{2TlhY 

(6-13b) 
They show clearly that the gauge-invariant Wigner representations of — i?iVi — ^ A(l) 

and ih-M 6^4(1) are given respectively by p and e. 

Using Eqs. ()6-8p and ()6-13p . we can express the Wigner representation of Eq. 
(|3-49p in a concise form as 

,^2 



2m 



K/ 



+ /X l-i:^(pe,rt 



=»K 



G^(pe,rt) = l, 



(6-14) 
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where the Groenewold-Moyal product is defined now by Eq. (|6-9p . 

From this point, we can directly follow the procedure described in ^4.31 Partic- 
ularly for the weak interaction case, we can adopt the transformation shown in ^5.11 
to an excellent approximation. We thereby obtain the Boltzmann equation in static 
electromagnetic fields as 

^^ +v-^ + e(E + -vxB] .^ = I^[f], (6-15) 



dt dr \ c J dp 

where v = p/m, and / is defined by Eq. ()5-3p . It is worth noting once again 
that, for discussing the electronic conductivity, we must incorporate the impurity 
scattering or electron-phonon interactions besides the two-body interaction in the 
collision integral Ip[f]. 

We finally comment on an alternative method that has been used frequently 
for deriving gauge-invariant transport equations .'^''^''^ It consists of adopting the 
ordinary Wigner transformation (j4-lap followed by the change of variables: 

p — >7v = p A, e — ;> 7r4 = e — 6^4, (6-16) 

c 

in the integrand. In the case of static electromagnetic fields, the corresponding 
Groenewold-Moyal product has been shown to be identical with Eq. ()6-9p .'^ There 
seem at least two difficulties in this change-of- variables method. First, it does not 
retain gauge invariance for electromagnetic fields with space-time variations. This 
may be realized by looking at the corresponding Wigner representation: 

G{p£, ri2ii2) = / (ffi2dh2 G(l, 2) e*i[P-(^/'=)^]*'i2-(^-^^4)*"}/^ (6-17) 

If we compare this expression with Eq. ()6-3ap . we realize that the phase /(1, 2) in the 
latter is replaced here by {e/c)A ■ fi2 — eA/^ii2. However, this transformation breaks 
gauge invariance through Eq. (|6-2p . except those cases where x(l) has linear depen- 
dence in (ri,ti). Moreover, it is unclear which space-time point we should adopt 
as the argument of electromagnetic potentials in Eq. (j6-17p . An exceptional case is 
when electromagnetic fields are static, where Eq. (|6T7p yields the same result as Eq. 
()6-3p . This difficulty was already pointed out by Stratonovich^^ ' and also discussed 
in detail by Serimaa et al.—'' Second, the change-of-variables method is apparently 
not applicable to superconductors where we have new types of functions, such as the 
pair potential Z\(l,2) and the anomalous Green's function F{1,2).^^'' The pair po- 
tential A{1,2), for example, transforms as ^(1,2) — )• exp{i^[x(l) + x(2)]} ^(1, 2) 
for Eq. (|6-2p . thus behaving like a wave function of charge 2e for 2 = 1. If one han- 
dles the product Z\(l,3)-F*(3, 2) using the change-of-variables method, one will end 
up having the operators d/dr and d/dt in front of A(p£, rt), thus failing to describe 
the fact that A{p£,rt) behaves like an effective wave function of charge 2e. 

In contrast, the transformation of Eq. (|6-3p is free from such difficulties and 
can be extended straightforwardly to superconductors. —' For the pair potential, for 
example, we only need to make use of the phase factor exp[— il(ri,^) — iI{f2,R)] 
[R = ^(ri+r2)] instead of e~*^(^'^' in Eq. ()6-3bp . The resultant Wigner representation 
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A{pe,rt) can be shown to transform as A{pe,rt) — t- exp[i^x(rt)]Z\(pe,rt) for 
Eq. (j6-2p . Moreover, the differential operators d/dr — i(2e/hc)A{rt) and d/dt + 
i[2e/h)Aii{rt) naturally appear in front of A[pe,rt). 

§7. Two-body correlations 

Here, we show that, with the self-consistent approximation considered in 5331 
we can also calculate two-body and higher-order correlations for a given <P. Thus, 
choosing a definite ^ amounts to determining the whole BBGKY hierarchy^ in the 
^-derivable approximation. 

To be specific, we consider the following two-body correlation: 

jC^m^l2,M) ^ (-0 (rcV^H(r)V;H(3'=)^li(4')V^Ii(2^)> - G*^(1,2)G'='(3,4). (7-1) 

The equation for /C is obtained by introducing an extra nonlocal potential Vl^(l , 2 ) 
on the Keldysh contour .'^''^ Let us regard 'Hit.) in Eq. ()3Tp as the unperturbed 
Hamiltonian. Then, the argument of ^2.21 tells us that all the effects of VF(1 , 2 ) 
can be incorporated into the S-matrix: 



•S'c = Tc exp 



2 

I 



Y,{-iy^' J dl I d24;l{l^)W'^{l,2)iPn{2n 



h 



(7-2) 



where the factor (—!)*"''•' stems from the transformation of Eq. ()3T3p . and the sub- 
script H signifies the representation defined by Eq. (|2-8p . Looking back at the second 
expression of Eq. ()3-lip . we also realize that Green's function with the extra potential 
W^^1,2)=W{1\ 2J) can be defined by 

Now, it is easy to see that correlation function (j7-ip is obtained from this Green's 
function as 






(7-4) 
w=d 



Here, the factor it is due to the permutation between ■011(2-') and 'i/'jj(2'-' ), and 
(—!)•'"*"•' cancels the factor (— 1)*+-' in Eq. (|7-2p . Equation ()7-4p may also be written 
concisely as IC = ±6G/6{f3Wf3). 

The connection (j7-4p enables us to derive the equation for IC, i.e., the Bethe- 
Salpeter equation, -' from Dyson's equation for G as follows. First, Dyson's equation 
is obtained from Eq. (j3-20p by adding an extra contribution of W. It may be written 
symbolically as 

G-^G = i, (7-5a) 
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G-'=(ih^^-Ky3-UW + IJ)f3. (7-5b) 

Hence, it follows that the first-order change 6G satisfies G^^SG + 5G^^G = 6, or 
equivalently, 

5G = G{-5G^^)G. (7-6) 

We also obtain 

6G S{t3Wt3) 

from Eq. ()7-5bp . Let us substitute this —6G~^ into Eq. (j7-6p . divide the resultant 
equation by S^f^Wfs), and set W = 0. We then express it in terms of Eq. (j7-4|) and 
the vertex: 

^ ' ^ ^h^ ' 6Gn' {2,2') h5&i'{l,l')5Gn'{2,2')' ^ ' 

where Eq. (j3-3ip has been used in the second equality. We thereby obtain the Bethe- 
Salpeter equation as 

V'' ^"' {ll\22') = ±G''J' {l,2')G^'' {2,1') ±ihY^ f d3 f d3' I dA f dA' G'''' {1,3') 

kk'W J J J J 

xG'=^'(3,l')r'=''=''''(3'3,4'4)/C"''-''^'(44',22'), (7-8) 



which determines IC for given G and F. 

We now reahze from Eqs. ([3^ . (fMTD . ([7^ . and dZS]) that the functional 
(p = ^[G] enables us to calculate one-particle, two-particle, and also higher-order 
correlations in a unified way. To put it another way, ^ specifies the approxima- 
tion definitely from one-particle to higher-order levels. This simple structure of the 
approximation is among the advantages of the self-consistent ^-derivable approxi- 
mation. 

The integral equation (|7-8|) can be solved formally as 

t = ±{l^ihGGry^GG, (7-9) 

where 1 and GG are defined by 

{lyi'^n' (11'^ 22') ^ 6.,j6,j>6{l, 2)6{l',2'), (7-lOa) 

{GGy'''^^'{n',22') = G'^'{l,2')ff''{2, l'). (7-lOb) 

It follows from Eq. ()7-7p that the vertex F has the symmetry: 

r'''"'{n', 22') = r^-'''**'(22', ll'). (7-lla) 

It also satisfies 

r^'>n'{n'^22') = -Yl {n)^k{fl)k'Afl)dn)l'A^'''''\^'h2'2)r, (7-llb) 

kk'W 
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which results from Eq. ()3-18p . These symmetries are useful in solving the Bethe- 
Salpeter equation. 

In H4:A\ we have adopted the local approximation for the self-energy. This local 
approximation may also be useful when calculating two-body correlations. Hence, 
we write down the Bethe-Salpeter equation in the local approximation. First, the 
vertex function is expanded as 

y^Qi[{Pi--r[-ei-t[)-{pi+-ri-ei+ti)+{p.2+-r'2-e2+t2)-{p2--r2-e2~t2)]/h 

(7-12) 

where pj± = pj it q/2, £j± = £j it uj/2, and rt is some space-time point around 
1,1', 2, 2'. Let us substitute Eqs. ()4Tb]l and (|7T2D into Eq. (US]). It then turns out 
that JC in the local approximation can also be expanded as 

K, ^^^^^^>- J ^27ThY J {27rhr J {27Thy'^ [PieuP2e2,quJ,rt) 

y^^i[{pi--ri-ei-ti)-{pi+-r[-£i+t[)+{p2+-r2-£2+t2)-{p2--r'2-£2-t2)]/!i^ 

(7-13) 

The equation for the Fourier coefficient /C" '-^-^ {piei,p2S2',(l^,ft) is formally the 
same as Eq. (|7-9p with the modifications: (i) Every matrix such as IC has the indices 
{pi£i,P2£2) to distinguish rows and columns, (ii) Integrations over internal variables 
are now given by J d^pd£/{2iTh)'^. (iii) The quantities 1 and GG are defined as 

{lf^"'{p£,p'£') = {2^h)H{p - p')6{£ - £')5ijSi,f, (7-14a) 

{GGf''^^\p£,p'£'; quo, rt) = {2T:hf8{p - p')5{£ - £')G'^' {p-£-,rt)&'' {p+£+, rt). 

(7- 14b) 

§8. ^-derivable approximation and conservation laws 

Careful consideration of conservation laws is necessary when tackling dynamical 
problems. For example, if an approximation you have adopted for a closed sys- 
tem does not conserve particle number, you will end up obtaining results that are 
completely nonsense. Thus, of crucial importance here is to find a general criterion 
with which conservation laws are satisfied. It is Baym^^- who presented a sufficient 
condition on this fundamental issue. To be more specific, Baym proved that con- 
servation laws are obeyed automatically in the (^-derivable approximation described 
in ^3.41 which still seems to be the only systematic approximation scheme with the 
essential property. In this section, we will provide a detailed proof of it, modifying 
Baym's original one with the equilibrium Matsubara formalisrrP^ onto the Keldysh 
contour .SSJ 
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8.1. Identities 

As preliminaries, we first derive various identities that are obeyed within the 
^-derivable approximation. Since this approximation becomes exact when ah the 
terms in Eq. (|3-3U|) are incorporated, those identities are also satisfied in the exact 
theory. 

First, we consider the fohowing gauge transformation: 



(5(2,1) 



eix(2)(5(2,l)e 



-ixW 



x(i) 



X(l) 




•1) 



Since it is composed of closed particle lines, the functional <P of Eq. (|3-3U|) is invariant 
through Eq. ()8-ip . including any approximate one with a partial summation. Thus, 
there is no change in ^ also at first order in x as 5^ = 0. By noting that ^ is a 
functional of G with the property ()3-3ip . the condition S<P = reads 



^2) 



dl / (i2Trf3i;(l,2)f3<5G(2,l) = 0. 



On the other hand, 6G can be written down explicitly from Eq. ()8-ip as 

5G(2, 1) = z [x(2)G'(2, 1) - G(2, l)x(l)] • (8-3) 

Substituting Eq. ()8-3p into Eq. ()8-2p and noting that x(l) is arbitrary, we obtain 

d2Tr^^[fsIJ{l,2)fsG{2,l) - G{l,2)fsIJ{2,l)f3] =0. (8-4a) 

It is further transformed into 

/ d2 [U^{1, 2)G^\2, 1) + U^^{1, 2)G^(2, 1) 

-G^(l, 2)i;^2(2, 1) - G^\l, 2)i:^(2, 1)] = 0, (8-4b) 

where we have used Eqs. (IHTHl . ([^T^ . and (IH^ to express S^'^G'^^ - S^'^G'^^ = 
S G^^ _l_ ^12(^A^ Equation ()8-4bp is what results from the gauge degree of freedom. 
Second, we consider the following Galilean transformation: 



G(2,l) 



exp 



R{t2)-^2\ G(2,l)expr^r%i)l , R{t) 



R{t) 




(8-5) 

Since the transformation only shifts the boundary of the system by R{t), <P is in- 
variant through Eq. ()8-5p . Hence, Eq. ()8-2p also holds here. The first-order change 
of G is now given by 

5G{2, 1) = R{t2)-V2G{2, 1) + %i)- ViG'(2, 1). (8-6) 

Substituting Eq. (|8-6p into Eq. (|8-2p and noting that R{t) is arbitrary, we obtain 

JQ{l)d\i = Q, (8-7) 
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with 

(5(1) = ^ih ^^~^^' I d2 [U^{1, 2)Gi2(2, 1') + i;i2(i^ 2)G^{2, l') 

-G^il, 2)i;i2(2, 1') _ ^12(1, 2)i;A(2, 1')] ,,^,, (8-8) 

where terms with derivatives of self-energies are due to partial integrations. 

Third, we consider a change of variables on Ci given by t — )• 9{t) = t + ^{t). 
Green's function is transformed accordingly as 



G(2,l) ^C/(t2)G(r2^2,ri0i)?7(ti), 



U{t) 



{do/dtYi^ 





where the factor (dO/dt)^'^ is introduced to cancel dt/d9 in Eq. (j3-23p upon the 
change of variables, thereby keeping ^ invariant through the transformation. Hence, 
it follows that Eq. ()8-2p also holds in this case. The corresponding 5G is given by 



5&\2,l) = \5ji 






+<5, 



'ji 






G^'^(2,l). (8-10) 



Substituting Eq. ()8-10p into Eq. ()8-2p and noting that (p{t) is arbitrary, we obtain 



dt^ 



where (7^int(ii)) and Qe{^) are defined by 



ih 



Qe(l)dVi, 



11 



(?iint(ti)> = ± J I d\^ I d2[u''{l,2)G''{2,l) + U'\1,2)G^{2,1) 



12) 



+G^(1, 2)i7i2(2, 1) + G^\l, 2)17^(2, 1)] , 
Q,(l) = T^n-^^ I d2 [^''{1', 2)G'^{2, 1) + U'\l' , 2)G^{2, 1) 

+G^(1, 2)i;i2(2, 1') + G'\l, 2)17^(2, l')] ,,^,. (8-13) 



Equation ()8Tip is what results from Eq. ()8-9p . 

The quantity (?^int(ii)) of Eq. (|8T2p is the interaction energy of the system. To 
see this, let us write down the (1,2) component of Eq. (|3-20p : 

lh-^-ki]G^^{l,2)- I [i;^(l,3)G^2(3,2) + i;i2(^^3)(^A(3^2)] (i3 = 0, (8-14) 

where we have used Eqs. (IHTJD . (fH^ . and (1^^ to express E^^G^"^ - U^^G"^^ = 
S G^^ _j_ j;i2(^A_ There is an alternative method of deriving the equation for G^^. 
It is based on 

i/i J-V'h(I) = Z^^(ti, -oo)[V'(ri),7i(ti)]^/(ti, -oo), 
oil 



42 Takafumi KiTA 

which stems from V'h(I) = ^^(^i, —oo)ip{ri)U{ti, — oo) and Eq. ()2-4p . By calculating 
the commutation relation with Eq. (j3-ip , this equation of motion is transformed into 



(ih-^^ - ki]i>n{l) - /^(l - 1')V'1i(1')V'h(1')V;h(1) = 0. (8-15) 

Let us multiply Eq. ()8-15p by ^{i/h)tp^{2) from the left, take the thermodynamic 
average with density matrix (j3-8p . and use Eq. ()3-15bp . Comparing the resultant 
equation with Eq. ()8-14p . we obtain the identity: 

i I dl' y(l-l')(^li(2)Vili(l')^H(l')^H(l)> = ±y / [^""(l^ 3)Gi2(3, 2) 

+i;^2(l,3)G^(3,2)]d3. (8-16) 

Setting 2 = 1 and performing an integration over ri in Eq. ()8-16p . we obtain an 
expression of the interaction energy in terms of the self-energies. We further take its 
complex conjugate to average the resultant expression and the original one. Noting 
Eqs. (j3-17ap and (j3-44p . we conclude that Eq. ()8-12p is indeed the interaction energy. 

8.2. Conservation laws 

We now show that conservation laws are automatically obeyed in the (^-derivable 
approximation. First, taking complex conjugate of Eq. ()8-14p and using Eqs. p-17ap 
and (|3-44p . we obtain 

ih^ - i^a) G^\i, 2) - / [G^(l, 3)i;i2(3, 2) + ^^^(l, 3)17^(3, 2)] d^ = 0. 

(8-17) 
We then subtract Eq. ([807]) from Eq. (fSTip with Eq. (j3^ in mind, set 2 = 1, and 
make use of Eq. ()8-4bp . We thereby arrive at the particle conservation law: 

^ + Vrj(l) = 0, (8-18) 

where n(l) and ji(l) are the particle density and the flux density, respectively, defined 
by 



n{l) = ±ihG^'^{l,l), (8-19a) 

(8-19b) 



2V1-V2 „i2/ 



2m 



2=1 



Second, let us operate ^ih{'Vi—'V2)/2m on Eqs. ()8-14p and ()8-17p . subtract the 
latter from the former, and set 2 = 1. It then turns out that the flux density j(l) 
satisfies 

^j(l) + -Vi0K(l) + ^Vi?7(l) = -Qil), (8-20) 

oti m m m 

where Q(l) is given by Eq. (j8-8p . and the tensor (1) is defined by 



e>f,{l) = ^'-X-{yu-y2i){Vi,-V2,)G'^{l,2%^y (8-21) 
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We further integrate Eq. ()8-20p over the whole space of the system and make use of 
Eq. ([8^ . We thereby obtain 



J- 1 3{l)d\, = - [!^V^U{l)d'n. (8-22) 

dti J J m 

Multiplying Eq. ()8-22p by m yields the total momentum conservation law. 

We finally focus on the energy conservation law. Let us operate =pi/i^ and 

^ih-^ on Eqs. ()8-14p and ()8-17p . respectively, add the two equations subsequently, 
and set 2 = 1. We then obtain 



where Qe{^) is given by Eq. ()8-13p . and £^{1) and ^^(l) are defined by 

ih^. 






2m 



2=1 



ih^ f d „ d , 
2m\dti dt2 



T— I— V2 + ^Vi G^^(l,2] 



•23) 

•24) 

•25) 



2=1 



The quantity £ {!) above signifies the kinetic energy density. We further integrate 
Eq. (j8^23p over the whole space of the system and make use of Eqs. (jS^lip and (j8^18p . 
We thereby arrive at the total energy conservation law as 



_d_ 



£^{l)d'n + {nint{ti)) 



Jil)-ViUil)d^n. 



•26) 



Equations (|8-18p . (j8-22p . and (|8^26p are the main results of this subsection. They 
clearly show that conservation laws are obeyed in the (^-derivable approximation. 

8.3. Sum rule for the spectral function 

The spectral function A defined by Eq. ()4-2p satisfies sum rule ()4-3p . which 
stands for the equal-time commutation relation of ip and lip^ . Although obvious 
intuitively, it may be worth showing that this property holds true in the dynamics 
of the (^-derivable approximation. 

Besides Eq. (|8^14p . we make use of the (2,1) component of Eq. (|3^2Up for this 
purpose, which reads as 

d 



ih 



-^-ki]G^\l,2)- I [i;^(l, 3)^21(3,2) + i;2i(l,3)G^(3,2)] (i3 = 0, (8^27) 

where we have used Eqs. (IHTJD . (f345p . and (I3^i3]l to express 1721^11 _ jj'i'^G'^'^ = 
r^G^i + 1721 ^A Let us subtract Eq. (fSTlD from Eq. ([8^ and set ^2 = ti subse- 
quently. It then follows that the self-energy terms cancel in the resultant equation, 
which can be expressed with Eqs. ()3^3p and ()4^2p as 



oti 2m 



Ail,2] 



0. 



(8^28) 



i2=ti 
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Let us take complex conjugate of Eq. ()8-28p with A*{1, 2) = A{2, 1) in mind, subse- 
quently exchange arguments as ri o r2, and subtract the resultant equation from 
Eq. (|8-28|) . We thereby obtain an equation of motion for A{riti,r2ti) as 

—^ + ^ — -^^A{riti,r2ti) ^ '-— — ^ -A{riti,r2ti) = 0. 

(8-29) 
Now, suppose A{riti,r2ti) = 5{ri — r2) at ti. Then, one may see easily that the 
second and third terms in Eq. (j8-29p vanish, so that dA{riti,r2ti) / dti = at ti. 
Hence, we conclude that sum rule ()4-3p holds at any t>ti. 

8.4. Conservation laws in the local approximation 

Equations (j8-18p . ()8-20p . and (j8-23p . which describe time evolutions of the par- 
ticle, momentum, and energy densities, respectively, are of fundamental importance. 
We express them here in the Wigner representation with the first-order gradient 
expansion so as to be compatible with the transport Eqs. (j4-17p and (j4-19p . 

Substituting Eqs. ()4-lbp and ()4-5ap into Eq. ()8-19p . we can express n{rt) and 
j{rt) in terms of A and (p. Those expressions are free from approximations and 
satisfy Eq. (|8-18p . To be more specific, the particle density n{rt) and the fluid 
velocity v{rt)=j{rt)/n{rt) are given by 

/ct T) ds 
p^/l(pE,rt),>(pE,r(), (8-30a) 

Equation (|8-18p is expressed in terms of these quantities as 

Qfi 

— + V(n^) = 0. (8-31) 

We next consider Eq. (|8-20p . By substituting Eqs. ()4-lbp and (|4-5ap into it, Eq. 
(|8-2ip is transformed into 

&fj{rt) = h J ^^^^A{pe,rt)(l){pe,rt) = mn{rt)vi{rt)vj{rt) + n^{rt), 

(8-32) 
where nf^{rt) signifies the kinetic-energy part of the momentum flux density tensor 
in the coordinate system moving with the local fluid velocity, given explicitly by 

n^{rt) = hj^^^A{pe,rt)(P{pe,rt), p = p-mv. (8-33) 

We next consider Q(l) in Eq. ()8-20p . Substituting Eq. ()8-16p and its complex con- 
jugate into Eq. ()8-8p . we obtain 

Q(i)=-/^M ^^g^~''^^ ^li(i)V^i,(iOV^H(iOV^H(i)), 
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with 1' = r[ti in this expression. We next express ('0jj(1)V'pj(1')V'h(1')V'h(1)) = 
P2{'r'i — f'i, ^2 ^ ;^i) ^^d expand ^2 ^ from ri up to the first order in the Taylor 
series. We thereby obtain 

Q(l) = -ViilP(l), (8-34) 

where nfj{^) is defined by 

77^(1) = -i|dV-r-^^(^Jj(l+)^Jj(l.)V;H(l-)V;H(l+)), (8-35) 

with 1± = (riibr/2,ti). Noting Eq. ()8-34p . we reahze that Uf- should be estimated 
with the local approximation in the first-order gradient expansion scheme. Hence, 
we express Ilf- with JC in Eq. (|7-1|) and substitute Eqs. (|3-5|) . ()4-lbp . and (|7-13p into 
the resultant expression. Carrying out a partial integration over f subsequently, we 
obtain 

P {ih)"^ f (fgdio f d^pde f d^p'de' f qiqj dVq 



^i^^'^^^- 2 J {2TTh)^J (2'KhY J (27r/i)4 y'i^'i+ q dq 

X [/C^2,i2 (pg^ p/g/. q^ ^^) ^ {2Tihf5{q)5{io)G^'^ {pe, rt)G^^ {p'e', rt)] . (8-36) 

Now, let us substitute j = nv, Eq. (18^32]) . and Eq. (f8^3il> into Eq. (l8^ and make 
use of Eq. (j8-3ip . We thereby arrive at 

dv 1 VU 

— + vVv + Vn + = 0, (8-37) 

ot mn m 

where Ilij = U^ + Uf-. By using Eqs. (|8-33p and (j8-36p . the tensor Uij can be 
written explicitly as 

d^pde pipj 



n,,{rt) = ^j^^ ^A{pe,rt)cl){pe,rt) 



{iKf f d^qdio f d^pde f d^p'de' f qtqj dVq 



2 J {2-KhY J {2Tihf J (27r/i)4 V^ -^ q dq 
x[K}'^^^'^{pe,p'e']quj,rt) + {2T:h)^5{q)5{u)G^'^{pe,rt)G^'^{p'e' ,rt)\,{^-m) 

where p = p — tuv, G^'^ is given by Eq. (|4-5ap . and IC is the solution of the Bethe- 
Salpeter equation (|7-9p in the local approximation. Note n_ is symmetric. 

We finally derive the differential energy conservation law. Instead of Eq. (j8-23p . 
it is convenient for this purpose to start from the expression of the energy density: 

^(i)^|^v;-Vi(^Ji(i')^H(i))|i,=i + ^y"dMF(n-rl)(^t^(i)^lj(i')^H(i')v;H(i)), 

(8-39) 
where t'l = ti is implied. Let us differentiate Eq. ()8-39p with respect to ti and 
remove the time derivatives of the field operators in the resultant expression with 
Eq. (|8-15p . We also carry out the first-order gradient expansion for the interaction 
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term by the procedure of deriving Eq. (|8-34p . We thereby obtain the differential 
energy conservation law: 

^^ + V • Mrt) = -j{rt) . VU{rt). (8-40) 

Here, je denotes the energy flux density defined in terms of the operator 

krt) = ^(V- V')V^li(r't)^H(rt)|^,^^, (8-41) 



by 

Urt) = ±^(V-V')V-V'Gi2(rt,r't) 



r'=r 



+ ^jdfVif){i,l{r^t)Jir+t) 



XTPuir^t)) - \ \ S-ff^^^- [(V;^(r_t)j(r+t)^H(r-t)) + {^^(r^t) 



A J df 

xj{r^t)i>H{r+t))]. (8-42) 

Noting Eq. (|8-40p . we reahze that Eqs. (|8-39p and (|8-42p should be estimated with 
the local approximation in the first-order gradient expansion scheme. It is carried 
out by the procedure of deriving Eqs. ()8-32p and ()8-36p . We thereby obtain 

£ = -mnv^ + £i, (8-43a) 

j'e = -mnv^v + £\v + IIv + Jq. (8-43b) 

Here, 77 is given by Eq. (j8-38p . and £\ and jq are defined in terms of p=p—mv by 

/d f) ds V 
J^^^,l^Mpe,rmpe,rt) 

±'-^ I ^^,[Ij''ipe,rt)G'^ipe,rt) + U'Hpe,rt)G^ipe,rt)], (8-44a) 

/d J) ds v 

{ihy r d^qdio r d^pds [ d?P'de' I ( q ■ P dVq 



+ 2 J (27rn)4 J (2^/i)4 j (27r/i)4 m ^^^^ + ^" 

X [K}'^^^'^{pe,p'e';quj,rt) + {2Tihf5{q)5{uj)G^'^{pe,rt)G^'^{p'e',r{^]^h) 



respectively. In deriving Eq. ()8-44ap . we have used Eq. (|8-16p . The quantity 8\ 
signifies the internal-energy density in the coordinate system moving with the local 
fluid velocity i), and jq denotes the heat flux density. 

Let us substitute Eq. ()8-43p into Eq. ()8-40p and transform the resultant ex- 
pression with Eqs. (|8-3ip and ()8-37p . We thereby obtain the differential energy 
conservation law as 

^ + V.(f,^+JQ) + 5;77,,|^ = 0, (8-45) 
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which for v = reduces to the heat equation. 

Equations (|8-3ip . (|8-37p . and (j8-45p are the differential particle, momentum, 
and energy conservation laws in the first-order gradient expansion, respectively. The 
quantities in those equations are defined by Eqs. (|8-30p . ()8-38p . and ()8-44p . They 
can be calculated microscopically with the solutions of Eqs. (j4-17p . ()4T9p . and (j7-9p 
in the local approximation. 

§9. Transport equation and fluid mechanics 

In most cases, fluid mechanics is constructed phenomenologically on the basis of 
the local conservation laws of particle number, momentum, and energj^^ without 
recourse to statistical mechanics. The Navier-Stokes equation thereby obtained is an 
archetype of nonlinear evolution equations embracing a wide range of nonequilibrium 
phenomena such as turbulence, chaos, and pattern formation. Since it consists of 
many particles, however, fluids should be described on a microscopic ground in a 
statistical mechanical manner, including its fluctuations. 

Taking a dilute classical monatomic gas as an example, we show here that the 
basic equations of fluid mechanics can be derived microscopically from the transport 
equation. To be specific, we start from the Boltzmann equation and subsequently 
adopt Enskog's expansion from local equilibrium^'-^ to derive equations of fluid 
mechanics appropriate for dilute classical gases. The consideration here enables us (i) 
to provide a microscopic basis of those equations, (ii) to calculate various transport 
coefficients such as the viscosity coefficient from first principles, and also (iii) to 
determine the form of the nonequilibrium distribution function. The distribution 
function thus obtained may be used with the method of ^ to calculate two-body 
correlations and fluctuations in nonequilibrium fluids. Thus, nonequilibrium systems 
are not mysterious at all, but can be handled in a statistical mechanical manner with 
the transport equation. One may also convince oneself after reading this section that 
typical nonequilibrium phenomena such as turbulence, chaos, and pattern formation 
in fluid mechanics lie fairly close to equilibrium. 

The contents of this section are described in detail in the classic textbooks by 
Chapman and Cowling^' and by Hirschfelder et al.®^' It may be useful, however, to 
include them here in the general framework of nonequilibrium quantum field theory, 
which can describe quantum effects and strong correlations far beyond the Boltzmann 
equation, to elucidate its great potential. Moreover, we proceed here by regarding 
Enskog's expansion as an approximate method to determine the distribution func- 
tion, i.e., the key quantity in statistical mechanics, rather than jumping directly onto 
calculations of transport coefficients.'^ The derivation here follows that in Ref. I83p . 

Our consideration starts from the Boltzmann Eq. (|5-4p . i.e., 

dt m dr dr dp ^ 

The collision integral Ip[f] in the dilute classical limit is obtained from Eq. (|5-4bp 
by neglecting the exchange effect as ±Vp-p^V*_p^ — )• and taking the classical limit 
as (1 lb /) — 7- 1. Performing a change of variables P2 = P + Q, P3 = P — q'/^j ^-nd 
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P4 = P + q' /2 subsequently, we obtain Ip[f] as 



/do Q f 
/2TTh)^ mj^^ [/p+(q-q')/2/p+(q+<7')/2 " fpfp+q]- (9-2) 



l2 



Here, da = dfiqi f dq'6{q' — q) [w'|V(q/_q)/2l/47r?i^] is the differential cross section 
in the coordinate system moving with the mass center^ and df2q' signifies the 
infinitesimal solid angle along q'. In the special case of the contact potential with 
Vq = const, the differential cross section can be expressed concisely in terms of the 
scattering length a = m\V\/A'irfi? by 

da = a^dflq, / dqb{(i - q). (9-3) 

The expression corresponds to the scattering of two rigid particles with radius a^ 
We will adopt Eq. (j9-3p below for the differential cross section. 

We next write down the conservation laws obtained from the Boltzmann equa- 
tion. Relevant quantities for this purpose are the particle density n{r,t), fluid ve- 
locity v(r,t), local temperature T(r,t), momentum flux density tensor n_{r,t) in 
the coordinate system moving with v{r,t), and heat flux density JQ{r,t). They are 
defined by 

d^p 

{2-Khf 



<r,t)= I j^:-^f{p,r,t), (9-4a) 



2 f (fp 



1 



with p = p — mv. These can also be obtained from general expressions ()8-30ap . 
(|8-30bp . (|8-44ap . (|8-38p . and ()8-44bp by neglecting the interaction, rewriting the 
resultant expressions with / in Eq. (|5-3p . and approximating the internal energy 
density by £i = ^nkBT for ideal monatomic gases. 

With these preliminaries, we now derive particle number, momentum, and en- 
ergy conservation laws. Let us multiply Eq. ()9-ip by 1, p/m, and p^ /2m, and carry 
out an integration over p. It then turns out that the contributions of collision inte- 
gral (|9-2p vanish in the resultant equations due to particle number, momentum, and 
energy conservations in the two-particle collision. Thus, the conservation laws are 
expressed with the quantities in Eq. (|9-4p as 



— + V{nv) = 0, (9-5a) 
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dv 1 VU 

— + V-VV + V77 + = 0, (9-5b) 

ot mn m 

-nk^ i-Q^ + v-VTj + V-JQ + n:Vv = 0, (9-5c) 

where A:B^ denotes the tensor product: 

A:B = Y^AijBji. (9-6) 

Equations ()9-5ap - ()9-5cp are the classical-gas versions of the general conservation laws 
(j8-3ip . (|8-37p . and (|8-45p in the local approximation, respectively. Indeed, Eqs. (|9-5ap 
and ()9-5bp are exactly identical with Eqs. ()8-3ip and ()8-37p . respectively, and Eq. 
(|9-5cp can be derived from Eq. ()8-45p by setting £[ = ^nk^T. 

On the basis of Enskog's expansion from the local equilibrium,'^''^ we now 
reduce the Boltzmann Eq. (|9-ip . which is nonlinear with the variables {p,r,t), to 
the problem of solving Eq. (j9-5p with the variables {r,t). Let us expand / as 

/(p,r,t) = /(i'=)(p,r,t)[l + ^«(p,r,t) + •••], (9-7) 

where /''^' denotes the local-equilibrium distribution: 

;<.e. = ,J?!I%,exp(-^V (9.8) 



(27rm/cBr)3/2 ^\^ 2mkBT 

with p = p — mv. This form of f^^^' has been chosen to satisfyP^''^ (i) the 
local-equilibrium condition that Ip[f] vanishes, and (ii) Eqs. (|9-4ap - ()9-4cp by itself. 
Hence, it follows that the higher-order corrections 99"' (j = 1,2,---) in Eq. (|9-7p 
should satisfy 

J^pnf(io)^U)=o, (n = 0,l,2) (9-9) 

There are five unknown functions in f^', i.e., n{r,t), v{r,t), and T{r,t), which are 
the same in number as Eq. (j9-5p . However, the latter still contains two unknown 
functions 77 and jq . 

We now express 77 and jq in terms of {n,v,T) so that Eq. (j9-5p becomes a 
closed equation for {n,v,T). Let us substitute Eq. (|9-7p into Eq. (|9-ip and regard 
the differential operators on the left-hand side as the first-order quantities in Enskog's 
expansion. Also noting 7p[/'^*')] = 0, we obtain the following first-order equation: 

^ + ii.^ + VC/.4^/<'«) =/»'[„<■)], (9.10) 

ot m or mk-Qi 

where I^p\ip'~^^] is given with Eqs. ([H^ and ([^^ by 

r(i)r (1)1 f '^^1 1 [^ A^e) Me) ( (1) , (1) (1) (1) 

(9-11 
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Using Eq. ()9-8p . we next express derivatives of f^ ' in Eq. ()9-10p with respect to those 
of {n, v,T). We further remove resultant time derivatives with the conservation laws 
in the local approximation. To this end, we make use of the following relations: 



77(le)=pi 



.■?'=". 



(9-12) 



with P = nk-QT and 1 as the unit tensor, which are obtained by substituting / = /''*') 
into Eq. ()9-4p . The left-hand side of Eq. ()9-10p is thereby transformed into 



9/(ic) p QfM p.^u 
m dr mk-QT 



dt 



+ 



(Ic) 



2 kk 



l]:Vv + [k 



5\ p 



2/771 



VlnT 



(9-13) 



where fc is a dimensionless vector defined by 



k =p/y/2mkBT. 
We now introduce a couple of additional dimensionless quantities: 



V 



2kBT 



-V, 



f = n ^'^r. 



and the mean-free path: 



m 



/ 



4\/27ra2n 



(9- 14a) 



(9-14b) 



(9-15) 



Let us rewrite Eqs. (j9-3p and (j9-8p in terms of these dimensionless quantities, put 
the resultant expressions into Eq. ()9-lip . and substitute it together with Eq. (|9-13p 
into Eq. ()9-10p . We thereby obtain 



2^2 fd?q fd^q'6{q' 



■spKlv}!'^ J iir J 4tt 



~k^-{k+qf r fl) 



(1) ,J1) 



\}fk+{q-q')/2 + ^k+{q+q')/2 "^k ^k+ 



2 kk- 



k^ 



l\:Vv + \k 



fc-VlnT 



(9-16) 



where V = d/dr, and we have redefined ^p^^' as a function of fc = p/\j2mk-BT. 
Similarly, Eq. (j9-9p is transformed into 



/d3A;e-'='fc>^^^ = 0. (n = 0,l,2) 



(9-17) 



Equation ()9-16p may be regarded as an integral equation for 99^ with condition 

The right-hand side of Eq. ()9-16p suggests that we may solve the equation by 
adopting the expression J^'l^ 



^'^ 



-In'/' 



A'^/^k^) (kk - ^l\ :Vv + A'/\e)k-VlnT 



(9-18) 
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Table I. Definitions of tensors in Eq. H9-19|) . wiiere 5'o (e) ~ 1 and Si{e) = 1 + a — e are Sonine's 
polynomials!^''^ 



Q 


w^ 


ni 


Tfe 


5/2 
3/2 


fefc-(fcV3)l 
k 




A^k^)yvZ 



where ^^'^ and ^'"^ are two unknown functions; they are distinguished here by the 
superscripts a = 5/2, 3/2 for the reason that will be clear shortly. Substituting Eq. 
(|9-18|) into Eq. (|9-16p . we obtain a couple of separate equations for A^ as 

2^2 [ d^q f d^q' 5{(/-ql ^^2_^^^ )2 r 

V^ y 17 7 17 q [ Jfe + ^fe+q - ^k+{<l-'l')/2 - ^k+{q+q')/2 

= e-^^nl, (9-19) 

where 7^^ and T^ are tensors given in Table [H The factor e~^ on the right-hand side 
is identical as a function of k"^ with the weight function of Sonine's polynomial.^^ ''S^ 
Hence, we solve Eq. ()9-19p by expanding ^"(e) in terms of Sonine's polynomials 
{S1{e)} as 

oo 

A^{e) = Y,o^^S1{e). (9-20) 

£=0 

By using the orthogonality of Sf{e) as well as Eqs. ()9-18p and ()9-20p . constraint 
(j9-17p is then transformed into 

al''^ = 0. (9-21) 

Hence, we will remove the ^ = term in expansion (|9-2Up for a = 3/2 hereafter. Next, 
we take the tensor (a = 5/2) or the vector (a = 3/2) product of Eq. ()9T9p with 
Sf {k'^)yV^ / A-K and subsequently perform the integration over k. Equation (j9-19p is 
thereby transformed into a set of linear algebraic equations for a^ as 

5^^,",,a?,=7^^ (9-22) 

Here, 7^" is given by 
d^k 



n?^J^e-'^'snk'){w^,n) 



5 

-V7r6eo -.a = 5/2 

15 ' ^^-^^^ 

-—\/Tr5ei : a = 3/2 
lb 



with (>V^,7^^) denoting W^ : 111 and W^ • Ul for a = 5/2 and 3/2, respectively. 
Besides, TZ, can be expressed with a change of variables k ^ k — q/2 as 



1 ff/ — 



/•OO /"OO /"(X) 

^^ / dke-^^'k" dqe-'^'l^q^ dq'6{q'-q) 
x[l^,,{k,q,q) + ire>ik,q,-q)-2I^,,{k,q,q')], (9-24) 
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Table II. Values of a? 


obtained by solving Eq. 


(|9-22|l. 


(y. dp 0-1 


a a 
"2 "3 


af 


5/2 2.2511 0.1390 
3/2 -1.7036 


0.0233 0.0058 
-0.1626 -0.0371 


0.0018 
-0.0117 



where I^pik, q,q') is defined by 

(9-25) 
and I^, {k, q, ±q) is obtained from Eq. (|9-25p by tlie replacement k — q' /2 — ;• fc =F q/2- 

It is easy to calculate the first few elements of Eq. ()9-24p analytically; we obtain 
Tq/^ ^ j,3/2 ^ ^^ ^5/2 ^ ^3^/2 ^ _^^^^ ^5/2 ^ 205/48, and T^/^ = 45/16. Besides, 
those for a general £^' may be evaluated numerically without difficulty. Now, we can 
solve Eq. (|9-22|) approximately using T^^ and TZ^ of ^'^' ^ ^ci and the errors can 
be estimated by increasing the cutoff £c- Table |TI] shows the values of a" thereby 
obtained. The convergence in terms of ic turns out to be quite rapid: for example, 
the values of Oq and a^ listed in Table |lT] are different from their lowest-order 
values 5-y/7r/4 (^c = 0) and —15^/^T/lQ (£c = l), respectively, by only ~ 2%. Thus, the 
expansion in Sonine's polynomial is quite efficient in obtaining the solution. 

Let us substitute Eqs. ([921), ([931), and ([948]) into Eqs. (f94dll and ([9^ . We 
then obtain expressions of the momentum flux density tensor and the heat flux 
density as 

dT 

JQ = -nCpK—. (9-26b) 

Here, P = nk^T denotes the pressure for the monatomic gas, Cp = |/cb is the specific 
heat per particle under constant pressure, and v and k are the kinematic viscosity 
and thermal diffusivity, respectively. The latter two quantities are defined in terms 
of the mean-free path / in Eq. (j9-15p as well as Og and a^ in Table |TI] as 



1/ /2fcBT^5/2 

4 V m 



l\l^—aQ, (9-27a) 



which clearly have the same dimension. Their ratio Pr = v j k constitutes the Prandtl 
number,'^ i.e., one of the fundamental dimensionless quantities of fluid mechanics. 
Our theoretical estimate Pr = 0.66 from Table HI] and Eq. (|9-27p agrees excellently 
with the experimental value 0.67 of Ar at T = 273 K and P = \ atm, and also those 
of many monatomic gases. "^ 

Thus, we have succeeded in expressing 77 and jg in terms of {n,v,T) as Eq. 
(|9-26p . With these expressions, Eq. (|9-5p now forms a closed set of equations for 
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{n,v,T). Among them, Eq. ()9-5bp is the standard Navier-Stokes equationp^ and 
Eq. (|9-5cp with t? = is the heat equation. It is worth emphasizing once again that 
these equations have been obtained by carrying out Enskog's expansion from local 
equilibrium up to the first order. Hence, it follows that, within their applicable range, 
the deviation of any solution to Eqs. (|9-5p and (j9-26p from the local equilibrium 
should be small to be expressed as Eq. dMI with Eqs. dEH]), (|9a8D . (|?^20D . and 
Table |lll Noting that those nonequilibrium phenomena are well described using 
the equations of fluid mechanics, one may realize that most of the energy in the 
turbulence, pattern formation, etc., is still dominated by heat, which can be described 
appropriately by equilibrium statistical mechanics. We finally point out once again 
that, using Eq. (|9-7p thereby obtained, we can calculate fluctuations in the fluid with 
the method described in ^ 

§10. Summary and future problems 

We have described the fundamentals of nonequilibrium quantum field theory on 
the basis of a many-body perturbation expansion with the Keldysh Green's function 
in a self-contained manner. The time-evolution operator and S-matrix were intro- 
duced concisely in §2; using the latter, we can calculate an expectation value of an 
arbitrary operator perturbatively with Eq. (|2-19p . It was shown in §§3.1- 3.3 that 
the perturbation expansion can be carried out with slight modifications of the equi- 
librium Matsubara formalism, specifically in the integration contour and Feynman 
rules. The formulation given here for a nonrelativistic two-body interaction can be 
applied easily to other many-body systems by merely changing the starting Hamil- 
tonian so as to be suitable for the relevant problem. A special consideration was 
given in §3.4 to the self-consistent perturbation expansion, i.e., Baym's ^-derivable 
approximation, of satisfying various conservation laws. It consists of solving Eqs. 
()3-20p and ()3-3ip self-consistently for an appropriately chosen functional ^. These 
equations form a firm microscopic basis for tracing nonequilibrium time evolutions 
of various many-body systems. Indeed, active investigations in high-energy physics 
have been carried out over the last decade to solve those equations numerically with- 
out any further approximations.'^T^ 

On the other hand, coupled equations ()3-20p and ()3-3ip may be simplified fur- 
ther into first-order differential equations in the phase space, i.e., quantum transport 
equations (§§4.1- 4.4), by transforming them into Wigner representations and trun- 
cating the resultant gradient expansion at the first order. The procedure may be 
justified if the space-time inhomogeneity occurs over a length or time much larger 
than the microscopic length or time characteristic of the relevant system, e.g., average 
interparticle spacing. Indeed, this condition is met in a wide range of nonequilib- 
rium phenomena. A wide applicability of transport equations was also confirmed by 
a numerical simulation.'^ Thus, quantum transport equations ()4-17p and (|4T9p in 
the phase space may be used as a convenient alternative to Eqs. (|3-20p and (j3-3ip 
in the full coordinate space. The reduced equations in the phase space also enabled 
us to obtain an expression of nonequilibrium entropy as Eq. ()4-27p . which evolves 
as Eq. ()4-26p and also embraces the Boltzmann entropy for dilute classical gases 
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as a high-temperature weak-couphng hmit (§5.1). This may be regarded as a defi- 
nite contribution to the basic issue of how and in what circumstances we can define 
nonequihbrium entropy. We also showed that the expression satisfies the law of in- 
crease in entropy within the self-consistent second-order perturbation expansion; it 
remains to be clarified whether the proof can be extended beyond the second order. 
Concerning the law of increase in entropy, we presented an argument at the end of 
§3.1 on why thermalization is achieved in this formalism by only evolving the system 
mechanically. 

In §5, we performed further reductions of relevant variables with a couple of 
different approximation schemes, i.e., the quasiparticle and quasiclassical approxi- 
mations. The quasiparticle approximation was shown to reproduce the Boltzmann 
equation. Sections 6-8 were devoted to quantum transport equations for a charged 
system in electromagnetic fields, a derivation of the Bethe-Salpeter equation for 
two-body correlation functions in the ^-derivable approximation, and a proof of 
conservation laws in the ^-derivable approximation, respectively. Finally, in §9, we 
presented a derivation of the Navier-Stokes equation for a dilute classical gas from 
the Boltzmann equation described in §5.1. 

Thus, readers may have been convinced that the nonequihbrium quantum field 
theory provides us with a powerful theoretical tool for describing general nonequi- 
hbrium phenomena based on Hamiltonians. Particularly, the ^-derivable approxi- 
mation enables us to handle one-particle, two-particle, and higher-order correlation 
functions on an equal footing for a given 'P so as to satisfy conservation laws. To put 
it another way, choosing an appropriate ^ amounts to fixing the whole BBGKY hier- 
archy:—-* Moreover, higher-order correlations and collisions beyond the two-particle 
ones can be incorporated systematically in the self-consistent perturbation expan- 
sion. These clear and simple structures of the approximation scheme are among 
the definite advantages of the (^-derivable approximation. Since no such schemes 
seem available in classical statistical mechanics, it may be used even for describing 
classical phenomena as a starting point to take the high-temperature limit even- 
tually. Note also that the Boltzmann equation and Navier-Stokes equation can be 
derived from Eqs. ()3-20p and ()3-3ip with successive reductions of relevant variables. 
Thus, those phenomena that have been discussed traditionally using the Boltzmann 
or Navier-Stokes equation can be handled microscopically and systematically by 
the (^-derivable approximation of nonequihbrium quantum field theory. They in- 
clude two-body correlations of nonequihbrium fluids, especially near nonequihbrium 
phase transitions,'^''^ which are clearly beyond the scope of the Navier-Stokes 
equation.'^ Thus, typical nonequihbrium phenomena such as pattern formations, 
which have been discussed conventionally on the basis of deterministic evolution 
equations for averaged quantities (e.g., the Navier-Stokes equation for the average 
velocity) ,'I^'12li can also be treated using this approach so as to incorporate fluctu- 
ation effects with the Bethe-Salpeter equation. It also enables us to improve the 
Boltzmann equation systematically so as to include higher-order collisions and/or 
quantum effects. It is interesting to see how Enskog's expansion described in §9 may 
be extended to systems with strong correlations and/or quantum effects. 

The approach may also be applied to superfluid phases,!® especially to Bose- 
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Einstein condensates which have been actively investigatecP^T^ since the reahzation 
of Bose-Einstein condensation on an atomic gas in lOOS.!^ In this context, it is worth 
pointing out that an extension of the (p-derivable approximation to Bose-Einstein 
condensates has met a particular difiiculty^T^ of not being able to reproduce 
the gapless Nambu-Goldstone bosorl^T^ of broken U{1) symmetry, which is the 
Bogoliubov mode in the weak-coupling regime!^ It is only very recently that the 
task has been successfully carried out™ Thus, applications of this approach to 
Bose-Einstein condensates remain mostly unexplored. They include a microscopic 
derivation of the two-fluid equations^ with a definite temperature dependence of 
the viscosity coefficient, etc. See Ref. [87| ) for the present status of the issue. 

Having derived an expression of nonequilibrium entropy in 21 it is also in- 
teresting to seek its possible roles in nonequilibrium phenomena. Much effort has 
been directed towards extending equilibrium statistical mechanics to nonequilibrium 
situations .'^''^''^''^^ One of its central issues has been whether there are any 
universal probability distributions in nonequilibrium systems corresponding to the 
microcanonical, canonical, and grand canonical ensembles of equilibrium statistical 
mechanics, or to put it another way, whether there are any thermodynamic functions 
that take minimum in nonequilibrium steady states. Various proposals have been 
presented on this issue, but no definite answer seems to have been obtained yet. 
In this context, it may be worth pointing out that entropy is the key quantity in 
equilibrium statistical mechanics. Indeed, every equilibrium statistical ensemble can 
be identified as the maximum of the Gibbs entropy (jl-ip under some constraints, as 
shown by Jaynes.*^ Thus, it may be worth asking (i) how entropy changes as we 
drive the system out of equilibrium and (ii) how far into the nonequilibrium region 
the maximum property of entropy holds true in the space of steady states. It should 
be noted at the same time that investigations of nonequilibrium states may be car- 
ried out by solving some evolution equations, e.g., Eqs. (j3-2Up and (j3-3ip in this 
paper, without recourse to any thermodynamic function such as entropy. However, 
a thermodynamic function, if it does exist and takes minimum in nonequilibrium 

steady states, will be useful in clarifying and classifying nonequilibrium phase tran- 
sitions.lSJ'C© 

The above questions (i) and (ii) on entropy were investigated using the Rayleigh- 
Benard transition of a dilute classical monatomic gas.l^ It was shown that entropy 
as a function of mechanical variables becomes smaller under a heat fiow than without 
it, in accordance with equilibrium statistical mechanics. Further, entropy was shown 
to be larger with convection than that without it in the region where the convection 
is stabilized. Thus, entropy as a function of mechanical variables apparently takes 
its largest possible value in the example, implying that the principle of maximum 
entropy may hold even in steady states beyond equilibrium o More investigations 



*' It is worth pointing out that this consideration of nonequiHbrium entropy differs substantially 
from that of extended irreversible thermodynamics developed by Jou et al.'^^ in that: (i) we 
adopt an expression of nonequilibrium entropy with interaction; and (ii) we focus only on those 
steady states that are solutions to the Boltzmann equation (or more generally, Dyson's equation 
on the Keldysh contour). On the other hand, extended irreversible thermodynamics makes use 
of noninteracting entropy, adds variables characteristic of nonequilibrium to entropy, and directly 
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are required to establish the above conjecture, however. 

Acknowledgements 

The author would like to thank Dr. H. Yaniashita for a careful reading of the 
manuscript and resultant useful comments. 

Appendix A 

Second Quantization 

The second quantization method enables us to handle many-particle systems 
concisely and conveniently. It is completely equivalent to the description with many- 
body wave functions in the configuration space. However, this equivalence has not 
been explained adequately in the literature, thereby leaving students who want to 
learn many-body physics in an awkward situation. We explain here the above equiv- 
alence as concisely and clearly as possible starting from the symmetry in identical 
particles under permutation. This will be carried out without assuming that the A^- 
body wave function can be expressed as a superposition of products of N one-particle 
wave functions. 

A.l. Permutation 

We first explain the permutation in the group theory.l^lJ A permutation signifies 
an action that A^ persons sitting on A^ chairs in a circle exchange their positions. It 
can be expressed by the operator: 

P=( ' ' ' - ^■), (A.l) 

\ Pi P2 P3 ■■■ PN J 

where 1 < Pi < N with no duplication among pi's. Thus, the person previously 
on chair i has moved onto chair pi. The number of distinct permutations is easily 
identified to be A^!. 

Among the permutations, there is a special class called cyclic permutation. It 
corresponds to each person moving onto an adjacent chair systematically as 

C^dllZ 'l' ^)^(123....), (A.2) 

A cyclic permutation of 2 elements is called transposition, which is expressed as 

A2^( 2 I )=(12). (A.3) 

We now summarize theorems on the permutation. See, e.g., Ref. I107p for the 
proofs. First, every permutation can be expressed as a product of cyclic permuta- 



maximizes the entropy in terms of those variables without asking whether the solution satisfies the 
Boltzmann equation. Thus, states identified with their procedure may not correspond to real steady 
states in nature. 
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tions. For example, 



Pn 



12 3 4 5 6 
2 5 6 4 13 



(4) (3 6) (12 5) = (3 6) (12 5). 



(A-4) 



The identity permutation is often dropped as (4) above. Second, every cyclic per- 
mutation can be expressed as a product of transpositions as 



(12)(13) ••• (1A;-1)(1A:), 



where operations on the right-hand side proceed from right to left. Combining the 
two statements, we realize that every permutation can be expressed as a product of 
transpositions. For example, Eq. (|A-4p may be expressed as 



1 


2 . 


• k-2 


k-1 


k 


2 


3 •• 


■ k-1 


k 


1 



Pa = (3 6)(12)(15). 



Another example is given by 



12 3 4 5 6 
2 5 6 3 14 



(364) (125) = (36) (34) (12) (15). 



An odd (even) permutation is a permutation that can be expressed as a product of 
an odd (even) number of transpositions as Pa (A) above. It can be shown that every 
permutation is either odd or even. 

A. 2. Permutation symmetry in systems of identical particles 

Many-particle systems composed of identical particles have a special symmetry 
under permutations. Consider an A^-particle system described by the Hamiltonian: 



H 



N 

E 



2m 



+ U{r, 



N 






Ti 



E'?'+Ei: 



i(2) 

ij 1 



(A-5) 



i=i 



i<j 



where pj is the momentum operator, and U and V denote the one-body and two- 
body potentials, respectively. This Hamiltonian commutes with Eq. (jA-ip as 



PHP 



-1 



H. 



(A-6) 



which is exemplified for the case of A^ = 2 as follows: 



P12HP, 



-1 
12 



*2 *2 

T^ + T^ + U{r2) + [/(n) + V{\r2 - nj 
Zm Zm 



12^12 



P15P1 



H. 



Thus, any permutation upon H only changes the order of summation to keep H 
itself invariant. It follows from Eq. ()A-6p that P and H can be diagonalized simul- 
taneously.l^^ In cases where H has no time dependence, the expectation value of P 
does not change in time. 

To begin with, let us determine eigenvalues of transposition (|A-3p . Noting that 
P12 is the identity permutation, we conclude that an eigenvalue ep^^ of A2 satisfies 
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Ep =1. Hence, it follows that epjj can be either 1 or —1. To put it another way, a 
wave function of identical particles either remains invariant or changes its sign upon 
a transposition. It was Pauli in 1940, who identified the connection between the 
eigenvalue epjj and the spin, which is summarized as 

!■-./, \ f 0,1,2,- •• (integer) 

_ f +1 (boson) • • J 

^■^''~l -1 (fermion) ^^ ^'^'''- | i, -,-,.. • (half integer) 

For example, the electron, proton, and neutron are all fermions with spin 1/2, 
whereas the photon is a boson with spin 1. The rule applies also to composite 
particles. Thus, the hydrogen composed of one electron and one proton, whose total 
spin can be either or 1, behaves as a boson. 

The above consideration can be generalized easily to the case of general per- 
mutations. Remembering "every permutation can be expressed as a product of 



transpositions" and "every permutation is either odd or even" in ^A.ll we obtain 
the eigenvalue £p = (±1)^ of a permutation P as 

/,-.\P _ / 1 '■ even permutation , , 

1^ ±1 : odd permutation ' 

where the upper (lower) sign corresponds to bosons (fermions). 

We finally write down explicitly how P operates on an A^-body wave function in 
the configuration space. Let us express the wave function as 

^u{xi,X2,--- ,xn), (A-8) 

where v signifies a set of quantum numbers to specify the A-body wave function, 
and Xj = {rj,aj) (j = 1,2,- •• , N) with Vj and Oj denoting the space and spin 
coordinates, respectively. The operation of Eq. (jA-ip on I^tv is defined by 

P^u{xi,X2,--- ,xn) = ^u{xpi,Xp2,--- ,Xpj^). (A-9a) 

On the other hand, the fact of 'Pu being an eigenfunction of P with the eigenvalue 
(±1)^ is expressed as 

P'^u{xi,X2,--- ,xn) = {±lf'^u{xi,X2,--- ,xn). (A-9b) 

A. 3. Constructing the eigenspace of permutation 

Equation (|A-9p tells us that an eigenfunction of identical particles should be an 
eigenstate of every permutation P, i.e., it should be either symmetrized (bosons) or 
antisymmetrized (fermions). We construct here the eigenspace of P where many- 
body wave functions of identical particles belong to. 

Following the creation and annihilation operators of the harmonic oscillator ,1^^ 
we first introduce the operators tf^ and tJj', which satisfy the commutation relations: 

['4>{x),ijj\x')]^ = '4^{x)'4'\x') =F V'"^(x')^(x) = 6{x, x'), (A-lOa) 
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mx),i>{x')]^ = [i;Hx),i>Hx')]T = 0, (A-lOb) 

with 5{x,x') = 6{r — r')6aa'- We next define the brackets |0) and (0| through 

V'(x)|0)=0, = {tlj{x)\0))^ = {0\i;^{x), (0|0) = 1. (A-11) 

These are basic ingredients to construct the eigenspace of P. 

We next introduce the state vector |xiX2 • • • xn) and its Hermitian conjugate as 



XiX2 • • -XN) 



l=i^^{xi)ip^{x2) ■ ■ ■ i'HxNM, (A-12a) 



/m 



(xiX2 ■■■xn\ = |2;i2;2 • • • xn)^ = -y={0\ip{xN) ' ' ' '4){x2)^p{xi). (A-12b) 

The space spanned by Eq. (|A-12aP naturally forms the eigenspace of P. Indeed, using 

i^\x2)i^\xi) = ibV'^(xi)V'''"(a;2), one can show \x2Xix^ ■ ■ ■ x n) = '^\xiX2X^- ■ ■ xn) 
easily. The commutation relation may be used repeatedly to prove 

\xp^Xp^ ■ ■ ■ Xp^) = (±l)'^|xiX2 • • • xn)- (AT3) 

The vector |xiX2 • • • xn) satisfies the normalization condition: 

{x'ix'2 ■ ■ ■ x'j^,\xiX2 ■■■xn) = -J^ ^{±1)^ 6{x'i,Xp^)6{x'2,Xp^) ■ ■ ■ 5{x'n,Xp^). 

(A-M) 
To prove it, we make use of the identity: 

^p{x'^)^\xi)^\x2) ■ ■ ■ iP\xn) = (5(Xi, xi)^t(x2) • • • ^\xn) 

+ {±lf5{x\,X2)^\xi)^\x:i) ■ ■ ■ i>\xN) 

+ ••• 

+(±l)^^t(^^)^t(^2) • • • i>HxN)i^{x\), (A-15) 

to place the annihilation operator ip{x'-) to the left of |0) and subsequently use 
Tp{x'A\0) = 0. For example, Eq. (|AT4p in the case of A^ = A^' = 2 is shown as 

(x;x'2|xiX2) = ^^{0mx'2)lPix[)i^Hxi)lPHx2m 

= ^m{x'2)[S{x[,Xl)i^Hx2)±S{x[,X2)i^Hxi) 

+ {±l)^Hxi)7P^{x2)i^{x[)]\0) 
= ^[^{x'i,Xi)6{x'2,X2) ±(5(Xi,X2)(5(x2,Xi)]. 

Thus, we have constructed the eigenspace of P. 
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A. 4. Bracket vectors for many-body wave functions 

We now introduce the bracket vector for wave function (jA-SP with symmetry 
(lA^ by 

I'^u) = / dxi / dx2 ■■■ dxN |a:iX2 • • ■ xn)'I'v{xi,X2-,- ■ ■ ,xn), (A-16a) 

{'I'„\= dxi dx2--- j dxN {xiX2---XNWt{xi.,X2.,--- -.xn) = \^u)\ (A-16b) 



where the integration over Xj = {rj,aj) is composed of that over Vj and the sum- 
mation over Oj. The bracket vector has the properties: 



ip{xi)\^^) = VN dx2--- dx'^\x2---x'^)%{xi,X2,--- ,x'^), (A-17a) 



i^{x2)ip{xi)\^u) = y/N{N - 1) / dxg • • • / dx^ jxg • • • x'^)l^,,(xi,X2,X3, • • • ,x^), 

(A-17b) 

iP{xn) ■ ■ ■ Hx2)i'{xi)\1'u) = Vn^.\0)<F^{xi,X2,- ■ ■ ,xn)- (A-17c) 

Thus, an operation of ip{x) on I'P^) has the effect of extracting the argument x 
from I'l^u)- The proof of Eq. ()A-17p may be carried out straightforwardly using Eqs. 
(|A-9p and (JA-ISP together with changes of integration variables. It follows from Eqs. 
dMH), (|A-12bp . and (|A-17cl) that 

{xiX2---xn\'1'u) = '^u{xi,X2,--- ,xn). (A-18) 

Substituting Eq. (|A-18p into Eq. ()A-16ap . we obtain 

I'F^) = / dxi / dx2 ■■■ dxN |a;i2;2 • • ■xn){xiX2 ■ ■ -xnI^u), 
which is equivalent to 

/ dxi / dx2 ■ ■ ■ dx]\[ |xiX2 • • • X]\[){xiX2 ■ ■ ■ xnI = 1. (A-19) 

Thus, the vectors {|xiX2 • • • xat)} form a complete set for the eigenspace of P. 

Equation ()A-16p can also be used to symmetrize or antisymmetrize any wave 
function l^(xi, X2, • • • , xn) which is not an eigenstate of P. Indeed, we only need to 
construct 

\^) = An dx[ dx2--- / dx']^\x[x2 ■ ■ ■ x']^)§{x[,X2, • • • , x'^), (A-20) 

with An denoting some normalization constant. Then, it follows from Eq. (|A-13p 
that integrations over x'l, • • • , x^ in Eq. ()A-20p with \x'iX2- • • x'^) in the integrand 
extracts from 'P(x'i,X2, ■ ■ ■ ,x'j^) only the symmetric or antisymmetric contribution. 
Equation ()A-20p may be regarded as projecting lf'(xi, X2, • • • , xat) onto the eigenspace 
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of P. The wave function of Eq. ()A-20p in the configuration space is obtained with 
Eqs. (|X44ll and (|A48]l as 

^{xi,X2,--- ,Xn) = {xi,X2,--- ,XAr|a/) = —^{±l)^^{Xp^,Xp2,--- ,Xpj^), 

(A-21) 
which indeed satisfies Eq. (|A-9p . 

A. 5. Orthonormality and completeness of bracket vectors 

Let us assume that the wave function ^i,{xi,X2-, • • • , xn) with symmetry (jA-9p 
satisfies the orthonormahty and completeness given by 

j dxi dx2--- / dxN^*,{xi,X2,--- ,xn)^u{xi,X2,--- ,xn) = K'ui (A-22a) 

^!p'i.(xi,X2,--- ,XN)Ki.x'i,x'2r-- ,x'n) = —^{±l)^6{x[,Xp^)6{x2,Xp2)--- 
V ■ p 

x5{x'^,xp^), (A-22b) 

respectively. These relations can be expressed alternatively in terms of the bracket 
vectors in Eq. (JA-16P as 

{^,,\^,) = 5,,,, Yl \^-)i^'^\ = 1- (A-23) 

V 

The proof proceeds straightforwardly by substituting Eq. ()A-16p into Eq. ()A-23p and 
making use of Eqs. ^^KM . (|AT4)) . (IXT9|1 . and (lA^22]) . 

A. 6. Matrix elements of operators 

We now see that matrix elements of operators for identical particles can be 
expressed concisely in terms of \^^) and ^(x). 

Equation ()A-5p exemplifies that every one-operator H^^' and two-particle oper- 
ator H^"^' for systems of identical particles may be expressed generally as 

N 

HW^Y.^hf\ H^^Y.^. (A.24) 

The matrix elements of these operators between ^*, and ^y can be expressed equiv- 
alently with {^i,i\ and \I/u) as 



/ 



dxi--- / dxArtP'*,(xi,--- ,xn)H^^^^v{xi.,--- ,xn) 

dxi{<Fy>\i;Hxi)h^^^{xi)\^y), (A-25a) 

dxi ••• / dxN^^'ixi,--- ,XAr)ij(^)l^j.(xi,--- ,Xn) 

^ I dxi I dx2{^u'\i^\xi)ij\x2)hfi^{x2mxi)\^,). (A-25b) 
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The equalities can be shown to hold in an elementary way by substituting Eq. ()A-17p 
into the right-hand sides of Eqs. ()A-25aP and (|A-25bp and using Eqs. ()A-9P and 

A. 7. Second quantization as an equivalent alternative to first quantization 

The preceding consideration has clarified that there are a couple of equivalent 
descriptions for many-body systems composed of identical particles. 

The first one is called the method of first quantization that makes use of Hamil- 
tonian (jA-Sp and the wave function (jA-Sp in the configuration space. Wave functions 
have to be symmetrized or antisymmetrized as Eq. (|A-9p . They may be constructed 
so as to satisfy Eq. (|A-22p . 

An alternative description is called the second quantization method, where the 
Hamiltonian is expressed concisely as 



n 



dxii[}^{ 



Xi 



pI 



+ U{ri) 



i){xi) 



+- 



2m 

dxi I dx2'iij'^ {xij'ijj'^ {x2)V {\ri - r2\)i>{x2)'4'{xi). 



The corresponding state vector is given in terms of Eq. (|A-12p by 
\Wy) = dxi--- / dxN\xi---XN)'l^u{xi,--- ,xn), 



(A-26a) 



(A-26b) 



which satisfies Eq. (^^23]) . It follows from Eq. (^^25]) that Eq. (^^26]) of second 
quantization is equivalent to the description of first quantization. 

A. 8. Second quantization in the absence of interaction 

We now express the Hamiltonian and brackets of noninteracting systems in terms 
of one-particle eigenstates. They also form a convenient starting point for a pertur- 
bation expansion with respect to the interaction. 

Consider the Hamiltonian: 



iio = / i^H^i) 



2m 



+ Uiri) 



%l:{xi)dxi, 



(A-27) 



and assume that the following one-particle eigenvalue problem has been solved: 

Vk{xi) = ekifkixi), (A-28) 



2m 



where k specifies the one-particle eigenstate and e^ is the eigenvalue. We also assume 
that the eigenfunctions (pkix) = {x\k) satisfy the orthonormality relations: 



{k\k') = / (fl{xi)ifk'{xi)dxi = 5kk', 



(A-29a) 



^ipk{xi)ip*k{x2) = S{X1,X2). 



(A-29b) 
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We next expand the operator ijj{x) in terms of ipki^i) as 

ip{x) = ^Ck(pkix), ip^x) = ^cl(pl{x). (A-30) 

k k 

Using Eqs. (|A-10p and (|A-29p . one can show easily that c^ and c|. above obey the 
commutation relations: 

[ck,cl'h = ^kk', [ck,Ck']T = [4'4']t = 0- (A-31) 

Let us substitute Eq. (|X30l) into Eq. (|A^27l) and take Eqs. ([A^ and (|A-29aP into 
account. Then, Eq. ()A-27p is transformed into 

'Ho = Y.^kcick. (A-32) 

k 

Thus, we have obtained the diagonalized expression of ?^o- 

We next express bracket vectors for T-Lq in the same representation. To this 
end, we first derive an A^-particle eigenfunction of 'Hq in the configuration space. It 
may be constructed from a product of A^ one-particle wave functions Y\j=i{xj\kj) = 

'^u{xi,X2,--- ,xn) by Eq. (|A-2ip. i.e., 

K{xi,X2,--- ,xn) = —'^{±l)^{xi\kp^){x2\kp2)---{xN\kp^), (A-33) 

p 

where u denotes u = {ki, ^2, • • • , k^)- In the case of fermions with the factor (— 1) , 
the summation over permutations in Eq. ()A-33p is equivalently expressed as^^a 



1'P{xi,X2,--- ,xn) = ^det 



{xi\ki) ■■■ {xi\kN) 



{xN\ki) ■■■ {xN\kN) 



(A-34) 



which is known as the Slater determinant. It follows from the properties of the de- 

/■p\ 

terminant that l^j) vanishes when two of the columns or rows are identical. We 
thereby arrive at the Pauli exclusion principle that no two identical fermions can oc- 
cupy the same one-particle state or the same coordinate (including spin). Thus, the 
Pauli exclusion principle naturally results from the permutation symmetry of iden- 
tical particles. By noting that all the fcj's are different, the normalization condition 
is transformed as 

1= dx^... dxM'l^!f\x^r■■,x^)\' = ^^Y.T.(-^)'^'^''Ili^p'J''p.) 

p pi i=i 

/.(F)n2 at / .(F)n2 / .(F)n2 

- (Niy 1^2-,^ ' W'p'.v,- i^\\2 2^2J^ ^) ^p'p- jY! • 

p pi i=i p pi 



Hence, we obtain Ajy as 



^(J) = VM. (A-35) 
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For bosons, on the other hand, every state kj can accommodate multiple par- 
ticles. Consider specifically the case where there are rij particles in the state kj 
(i = l,2,--- ,£■ i<N),le., 

e 
1^ = (ki,--- ,ki,k2,--- ,k2,--- ,ke,--- ,ki), X^^j = ^- (A-36) 

The corresponding wave function is given by Eq. (]A-33P with the upper sign, where 
kp. denotes a permutation of u in Eq. (|A-36p . The normalization condition is trans- 
formed as 

,{B)n2 N 



P' p i=i 



/rf^i---/rf^iv|l^f)(xi,---,x^)P = -L^J]^J](fc,,|fcp; 






(iv!)2 -z^ii^ •'I p^' m 



-ni!n2! • • • n^!, 



where the third equality originates from multiplying the result for P' = (identity 
permutation) by A^!. Hence, we obtain the normalization constant as 



^ ^Jnl\n2\ ■ ■ -nil 



<^ = ^===^- (A-37) 



Now that we have obtained A^-body wave functions in the configuration space, 
we write down the corresponding bra- vector |l^v)- Consider first the case of fermions. 
Let us substitute Eq. ()A-35|1 into Eq. ()A-33|1 with the sign (—1)^, put the resultant 
expression into Eq. ()A-26bp . and transform it further with Eqs. (|A-29p - (|A-3ip . We 

(V) 

thereby obtain \^i) ) as 

I'^f^) = 44---4J0). (A.38a) 

This expression is much simpler than Eq. (|A-34p in the configuration space, and 
the sign change of the latter upon the exchange of columns ki o kj is algebraically 
realized in \^i, ) through the exchange of positions between c^ and c^. with Eq. 
(|A-3ip . Specifically, the A^-body ground state corresponds to the state composed of 
A^ lowest one-particle-energy states, which is sometimes called Fermi vacuum. 

As for the case of bosons, we substitute Eq. ()A-37p into Eq. ()A-33p with the 
sign (+1)^ and transform the resultant expression in the same way as the case of 
fermions above. We thereby obtain \^l ) corresponding to v in Eq. (|A-36p as 

|^(B)) = fhL^ . . . i^|o). (A.38b) 

\/ni\ \/ni\ 

This state is identical in expression to that of a harmonic oscillator with multiple 
frequencies.'^^ 
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It follows from Eq. (|A-3ip that c||0) = holds for fermions. Taking this fact 
into account, we can put Eqs. ()A-38aP and ()A-38bp into a single expression as 






\iiy,) = iJ^...iJM_\o), (A.39) 

\/ni'. \/n£'. 



One may now realize that the second quantization method enables us to describe 
many-particle systems with identical particles concisely and conveniently. 

A. 9. Density matrix in second quantization 

The density matrix is undoubtedly one of the most fundamental concepts in 
quantum statistical mechanics. We express here equilibrium density matrices in 
second quantization. 

Consider first the canonical ensemble. Its density matrix can be expressed in 
terms of Hamiltonian (|A-26ap . its eigenstate ()A-26bp . and the eigenvalue Ei, as 

p = J2e^(^-^¥u){^.l F^-^ln^e-^^% (A.40) 

u u 

where /? = l/k^T with k^ denoting the Boltzmann constant and T the temperature. 
By using Eq. ()A-40p . the expectation value of an Hermitian operator O is evaluated 
as 

{6} = Y,e^^^~^''H'^AOWv) = TrpO. (A-41) 

V 

For example, the expectation value of the density operator h(x) = tjj^ {x)tjj{x) is 
expressed as 

(n(x)) =^e^'-^-^^\^:,\7p^{x)7p{x)\^u) =TrpV'^(x)V^(x). (A-42) 

Next, the density matrix of the grand canonical ensemble is given by 
p = Y^e^^''-^+^^^\KN){KN\, f2 ^ -iln^e-'^^^"-'^^), (A-43) 

uN ^ uN 

where p is the chemical potential, and 

N = i;\x)tl^{x) dx = ^ clck (A-44) 

k 

denotes the particle-number operator. The second expression of Eq. (|A-44p is ob- 
tained by making use of Eqs. ()A-29aP and ()A-30p . We now focus specifically on 
the noninteracting case with Hamiltonian (|A-32p . Substituting Eqs. (|A-32p . (|A-39p . 
and ()A-44p into Eq. ()A-43p . we obtain the grand canonical density matrix without 
interaction as 



E»-'""-'""i")"("i ,^,„ 

^ _ ln\j, = — 



W = n " v^ «, ^ . \n)i,s^-jL.\0)t, (A.45) 



66 Takafumi KiTA 

where |0)fc denotes the vacuum of the one-particle state k defined by Ck\0)k = 0, and 
the n summation runs over n = 0, 1, 2, • • • (n = 0, 1) for bosons (fermions). 
Equation ()A-26aP tehs us that it may be useful to introduce 

p^^\xi,x[) = Tr p'4>^x[)ip{xi) = {ip''{x[)tl){xi)), (A-46a) 

p^^\xi,X2;x[,x'^) = {i>\x[)i>\x'^)i>{x2)^ixi)), (A-46b) 

which are called reduced density matrices. These expressions apply to both the 
canonical and grand canonical ensembles with either Eq. (|A-40p or (|A-43p used as 
the density matrix. The expectation value of every one-body (two-body) operator 
can be expressed with p^^' {p^"^')- For example, the statistical average of Eq. (]A-26ap 
can be written in terms of Eq. ()A-46p as 



{%) = / dxi 



2m 



^(1) 



{xi,x[] 



X-, =X\ 



+- / dxi I dx2V{\ri - r2\)p^'^\xi,X2;xi,X2). (A-47) 

It is worth pointing out that Eq. ()A-46p can also be expressed with Eqs. (|A-17p and 
(|A-40p or (|A-43p in terms of wave functions in the configuration space. In the case 
of the canonical ensemble, for example, they are given by 

p^^\xi,x[) = ^e^^^-^-'^N dx2--- dxNMxi,X2,--- ,XN)Kix'i,X2,- ■ ■ ,xn), 

v J J 

(A-48a) 
;,4) = J^e^^^-^'-^iVliV-l) j dx^--- fdxNMxi,X2,X3,--- ,xn) 

V J J 

xt^;(x;,x'2,X3,--- ,xiv). (A-48b) 

These reduced density matrices were used by Yang^^ to discuss the concept of 
off-diagonal long-range order in superfluids and superconductors. 

Appendix B 

Self- Consistent Perturbation Expansion 



P^'^\xi,X2;X'i,X2 



We explain here the self-consistent perturbation expansion in the equilibrium 
Matsubara formalism. Its advantages are summarized below Eq. (|3-32p . The con- 
tents here also form the basis to understand nonequilibrium quantum field theory 
developed in §32] and El 

B.l. Thermodynamic potential and Green's function 
We consider a system described by the Hamiltonian: 

n = no + n,nt, (b-i) 

where T-Lq and T-imt are defined by Eqs. ()3-2ap and ()3-2bp . respectively, with no time 
dependence in U of Eq. (|3-3p . Note that Eq. (j3-3p includes the chemical potential. We 
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drop here the spin degrees of freedom for simphcity once again. The thermodynamic 
potential is given bj^-ED'El 

f2 = -^lnTVe-'^^, (B-2) 

with /3 = I/ZcbT. To calculate f2 perturbatively in terms of Timt, let us introduce 
the scattering matrix S{(3) as 

5(/3) = e'^^^e-'^^. (B-3) 

It obeys the equation of motion: 

^^ = --Hint(/3)cS(/3), 5(0) = 1, (B-4) 

with 

?iint(/3)^e^^«?iinte-^^°. (B-5) 

Looking back at Eqs. ()2-3p and ()2-4ap . we can solve Eq. (|B-4p as 

oo p p 

5(/3) = l + V(-ir / dTn--- dTi'Hint{rn)---'Hint(.ri) 
„-i -^0 Jo 



n=l 

r /-/^ . 

T^exp - / dr^int(r) 
. Jo 



(B-6) 



where 7^ orders the Hamiltonians into the ascending sequence of the arguments 
from right to left with an extra sign change upon every permutation for fermion field 
operators!^''^ One may easily see that this expression with Eq. ()3-2bp is equivalent 
to Eq. (|3-25p . We next introduce the noninteracting Matsubara Green's function: 

where ip{l) = e^i^O'0(ri)e~^i'^o and tl^\l) = e^i^o^t(y.^)e~^i^" with 1 = rin, and 
(• • • )o denotes the average in terms of the density matrix pQ = e~^^°/Tre~'^^°. 
Note that "0^(1) is not a Hermitian conjugate of V'(l) in the equilibrium Matsubara 
formalism. '20 '122} Writing Eq. ()B-3p as e~^'^ = e~'^^''5(/3), we can transform Eq. 
(fB^ into 

n = no-^ ln(5(/3))o = ^0 - ^ [(5(/3))oc - l] , (B-8) 

where Hq is the noninteracting thermodynamic potential, and the subscript c denotes 
retaining only connected Feynman diagrams in the expansion of (5(/3))o with Gq. 
See, e.g., Ref. [31]) for the proof of the second equality that taking the logarithm of 
(5(/3))o amounts to retaining only connected diagrams of (5(/3))o — 1. 

We focus next on the one-particle properties. They can be evaluated with the 
Matsubara Green's function: 

G(l,2) = -(r.^H(l)V^Ji(2)>, (B-9a) 
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where V'h(I) = e'^i^Vi(ri)e^'^i^, '0fj(l) = e'^i^V'^(^i)e^^^^, and (• • • ) denotes the 
average in terms of the density matrix p = e~''^/Tre~^^. Indeed, Eq. (|B-9ap is a 
direct extension of the one-particle density matrix (|A-46aP to incorporate the r de- 
pendence. The Matsubara Green's function can also be expressed with 5(/3) aa^^'ESJ 

G(l, 2) = -(r.5(/3)^(l)V;t(2))oc, (B-9b) 

which can be calculated perturbatively with Gq!^''^ 

To derive an equation of motion for G, we differentiate V'h(I) = e'^^ 'ijj{ri)e~'^^ 
in terms of ri to obtain 

OTi 

We then substitute Eq. (|BTp with Eqs. ()3-2ap and ()3-2bp into the above equation 
and calculate the commutator. The result can be expressed as 



gV'H(l) 



-i^iV'H(l)- /'d2y(l-2)^|j(2)^H(2)V'H(l), (B-10) 

where Ki is defined by Eq. (|3-3p with no time dependence in U, and V{1 — 2) = 
6{ti — T2)V{ri — r2). Let us multiply Eq. (|B-10p by — ^Jj(l') from the right, op- 
erate 7^, average the resultant expression with p = e^" /Tve~^ , and make use 
of (7;%^V'h(1')) = ^<^-^hW^(iO) _ ^(;l^ ^/)_ Tj^g f^^^i ^gg^i^ p^^ ^g written in 
terms of Eq. ()B-9p as 

-^ - ^i) G'(i, 1') + y d2 y(i - 2)(r.^^(2)V;H(2)^H(l)^Ji(l')) = <^(1, !')• 

(B-ll) 
The equation of motion for Go is obtained from above by setting V = 0, which may 
be expressed alternatively as 



Go\i,1')={-^-kA6{1,1'). (B.12) 



On the other hand, inspection of the perturbation expansion of Eq. (|B-9bp shows 
that G(l,2) satisfies Dyson's equationPJ'ED'l33 

G(l, 1') = Go(l, 1') + f d2 f d2'Go(l, 2)17(2, 2')G(2', l'), (B-13a) 

where i7(2, 2') is the (proper) self-energy, which cannot be separated by cutting any 
single-particle line. The first few Feynman diagrams in the perturbation expansion 
for U are given in Fig.[8l where straight and wavy lines denote Gq and V, respectively. 
Let us multiply Eq. ()BT3aP by Gq (3, 1), carry out an integration over 1, and make 
use of Eq. (|B-12p . We thereby obtain an alternative expression of Eq. ()B-13aP as 

(-^ - K^ G(l, 1') - J 17(1, 2)G(2, 1') d2 = 5(1, l'). (B-13b) 
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+ A/vA + 



Fig. 8. First few Feynman diagrams for the proper self-energy. 



Comparing Eq. ()B-13bp with Eq. (jB-lip . we arrive at the identity: 

d2U{l,2)Gi2,l') = - [ d2V{l-2){Tri^l{2)i,ni2)4^Rmlil')). (B-14) 



We now set 1' = 1+ in Eq. ()B-14p . where the subscript + denotes the presence of an 
extra infinitesimal positive constant in ri to place the creation operator 'i/;pj(l') to 
the left of the annihilation operators. We then multiply the resultant equation by 
=pl/2/3 and perform an integration over 1. The result can be expressed in terms of 
Eq. (|3^ as 

1 



{-H 



int/ 



=F 



2/3 



dl / d2E{l,2)G{2,l')\^,^^^, 



(B-15) 



which will play an important role below. 



B.2. Renormalization of self-energy 

Let us focus on the last four diagrams of the second order in Fig. [HI We then 
notice that the parts enclosed by squares below have the effect of inserting the first- 
order self-energy into Gq to renormalize Gq ^ G {G is denoted by a thick line). 



+ 




6\/\/0 



The observation suggests an alternative way to obtain the self-energy: consider only 
skeleton diagrams, i.e., diagrams without self-energy insertions, in the expansion for 
the self-energy and put Gq — )• G. We call this expansion the renormalized skeleton 
diagram expansion for the self-energy, which can be expressed graphically as follows. 



c:^ + ^^ 



The corresponding analytic expression is given by 

i:(l, 1') = ^5(1, 1') j d2 V{1 - 2)G(2, 2) - y(l - l')G(l, l') 

+ f d2 f d2'V{l-2)V{l' -2')[±G{1,1')G{2,2')G{2',2) 
+G(1, 2')G'(2', 2)G(2, 1')] + • • • , (B-16) 
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where G(l, 1') for ri = r( means G(l, 1^). 

The renormahzed skeleton diagram expansion is a self-consistent perturbation 
expansion. Indeed, E here is given as a functional oi G as U = E[G], whereas G 
obeys Dyson's equation: 



+ 



so that it necessarily depends on U. Hence, it follows that G and U have to be 
determined self-consistently in the renormahzed skeleton diagram expansion. The 
first-order self-energy, denoted by U , is expressed graphically as follows. 

Using S^^ in Dyson's equation ()B-13p reproduces the self-consistent Hartree-Fock 
approximation, as first pointed out by Luttinger.^ He also suggested more general 
approximations of incorporating higher-order contributions in Eq. ()B-16P into the 
self-energy of Dyson's equation (jB-lSp .''^^^ 

B.3. Self-energy as a functional derivative of (P 

The renormahzed self-energy in Eq. ()B-16p can also be obtained as follows. Let 
us introduce the functional ^ = ^[G] as the contribution of the skeleton diagrams 
in the perturbation expansion for i? with the replacement Go — )■ G. It may be 
expressed graphically as follows. 

<D = o^xo + Q + Q Q + N/ + ... 

The corresponding analytic expression is obtained from Eq. (|B-8p as 

^ = —■o[('5(/3))oc — l]skeleton,Go^>G 

= ^ / '^l / dl' ^(1 - l')[G{l, 1)G(1', 1') ± G(l, 10G(1', 1)] 

~Jb t ^^ t '^^' /"^^ f d2'V{l-l')V{2-2')[G{l,2)G{2,l) 

xG(l', 2')G(2', 1') ± G(l, 2)G(2, l')G(l', 2')G(2', 1)] + • • • . (B-17) 

Comparing Eq. ()B-17p with Eq. ()B-16p . we notice that there is an additional G in 
each order with an extra factor — (ibl)/2n/3 in the nth-order terms. The factor l/2n 
reflects the fact that 2n G's are equivalent in the nth-order closed diagrams for ^, 
and ±1 originates from the fact that there is an extra closed particle line. With 
these observations, one may convince oneself order by order that the self-energy of 
Eq. ()B-16p is also obtained from <P above by 

This is an exact relation with ^ defined by Eq. ()BT7p . Equation (|B-18p also pro- 
vides us with a systematic self-consistent approximation scheme called "(^-derivable 
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approximation" by Baym.— ' It consists of (i) retaining only part of the series in 
Eq. ()B-17p and (ii) determining G and U self-consistently using Eqs. ()B-13P and 
(|B-18p . Although it is equivalent to considering some partial series in Eq. (|B-16p . 
the ^-derivable approximation has a clear advantage that giving <P also enables us to 
calculate two-particle and higher-order correlations, as shown in ^ It is also with 
^ that one can prove conservation laws concisely; see ^on this point. 

Appendix C 

Luttinger- Ward Functional 



The thermodynamic potential is defined by Eq. (|B-2p . which can also be ex- 
pressed as Eq. ()B-8p . Below, we provide it with a couple of alternative expressions. 
The latter expression, i.e., the Luttinger- Ward expression, ^^' is quite useful for clar- 
ifying its general properties as well as for carrying out practical calculations. 

Let us introduce a parameter A' into V. = T-Lq + T-Lmt of Eq. (|B-2p as Q\i = 
—I3~^ InTr e~^*^^"^'*' ^^nt). We next differentiate Q\i in terms of A' to obtain 



dQy Tie-^^y'Hir.t {X'n 



hit /A' 



d\' Tre-'^^A' A' 

where {• • •)\' denotes the grand canonical average with respect to Tix' = Ho + X'T^int- 
Let us integrate the above equation in terms of A' from to A with 'H\'=o = Tio in 
mind. We thereby obtain 

f2x = Qo + f ^^"^^'^'' dX'. (C-1) 

We finally substitute Eq. ()B-15p into Eq. (|C-ip to arrive at 

^A = ^oT^y ^ldljd2Eyil,2)Gyi2,U). (C-2) 

We need Sx and Gx over < A < 1 before using this expression. 

Luttinger and WarcF^ found an alternative expression for i? as a functional of 
the self-energy E. It is more convenient^ to rewrite their result as a functional of 
G as 

Q = ±^Tr[ln(-G^i + T) + EG] + <?. (C-3) 

Here, G denotes the matrix with the element G(l, 1'), the symbol Tr signifies 

TiA= f A{l,l+)dl, (C-4) 

and (p is defined by Eq. (]BT7p . This expression enables us to calculate f2 directly 
from G and U, and also to obtain some exact results on thermodynamic quantities. 
Expression ()C-3p has an important property of being stationary with respect to 
a variation of G as 

. 5G(1', 1)^-0 -' ^ 6G{V,l)-^-6G{V,l)_ 



6n 1 

±-Tr 



SG{V,1) (3 



+ 



5G{V,l) 
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= 0, (C-5) 

where we have made use of Eq. ()B-13bp . i.e. G~^ = G^^ - S, and Eq. (jBTS]) . Thus, 
i? is invariant through an infinitesimal variation in G. See Ref. H2|) by de Dominicis 
and Martin for a more general proof of the stationarity as well as its extensions to 
the superfluid phase. 

Using Eq. (jC-Sp . we can prove the equivalence of Eq. (jC-3|) with Eq. ()C-2p for 
A = 1. To this end, we express as f^x the quantity obtained from the right-hand side 
of Eq. (|C-3p by the replacement T^int — t- AT^int, i.e., 

nx ^ ±^Tr[ln(-Go-^ + E,) + r,Gx] + <Px. (C-6) 

There are two types of A dependences in ilx. The first one originates from the n 
interaction lines in the nth-order contribution to ^Ai and the other from the implicit 
dependence through Gx- Let us differentiate f2x with respect to A. Due to Eq. (|C-5p . 
however, we need not consider the latter dependence in the differentiation. Noting 
dX^/dX = nA"/A in the explicit dependence as well as Eqs. (|B-18p and (|C-2p . we 
obtain 

Thus, Hx and Hx satisfy the same first-order differential equation in terms of A. 

To complete our proof of f2x = ^x, we now show that the initial value Hq = 
ib/3^-'^Trln(— G(^ ) satisfies ^o = -^o- To this end, we adopt the representation of 
diagonalizing Eq. (|B-12p to expand it as 

^ oo 

£=— oo k 

where e^ denotes 2£-7t/(3 and {2i+ l)vr//3 for bosons and fermions, and ^k and (pkifi) 
are the eigenvalue and eigenfunction of the operator Ki = —(h?/2m)Vl + U{ri) — ^, 
respectively. Then, i7o is expressed as 



^0 = ±^ E E ""'°^ i^^^*^ - ^^^)- (c- 



"/? 



e k 



The factor e*^*''+ originates from the subscript + in Eq. ()C-4p . By noting that ie^ is 
a pole of the function 

/(-) = ^^' (C-9) 

with the residue ±/3~^, the summation over I can be transformed into an integration 
on the complex z plane as 

dz .^ e^°+ 
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Fig. 9. Contours on the complex z plane. 



where C denotes the contour in Fig. [9l With the factor e^^+ , we can further make 
use of Jordan's lemma to deform C into C" in Fig. O We then perform an integration 
by parts with (e^^=Fl)~"'^ = ^P~^'d£^'^{^-\-G~ )■ It then turns out that there is nuh 
contribution from the edges of C . Finally, using the residue theorem, we obtain 



^n 



±^5]ln(lTe-^«'=) = f2o. 



"/? 



(C-10) 



Equations ()C-7p and (jC-lOp enable us to conclude Q\ = Q\; setting A = 1, we arrive 
at Eq. dnH). 

Appendix D 

Expression of Equilibrium Entropy 



We derive here an exact expression of equilibrium entropy S = —dQ/dT from 
Eq. (|C>3D.B3-I53 

Let us expand Green's function and the self-energy as 



.. oo 

G(l,2) = - ^ G(ri,r2;ie,)e-^^^(^i--2), 



P 



e=-oo 

OD 



^(1,2) = i Yl ^iri,r2;iSi)e-''' 



(ti— T2) 



/3 



(D-la) 
(D-lb) 



e=-oo 



where ei = 2£7r//3 and (2£+l)7r//3 for bosons and fermions, respectively. Substituting 
these expressions into it and carrying out Tr over r, we can transform Eq. ()C-3p into 



^ = ±^Y1 '^^{^^ i— + — (^^^) ~ ^^^-] + — (^^^) G{iee)}e'^'°+ + <P. 

Here, X, G{ie£), and ^{i£i) signify matrices with the elements 

K{ri,r[) = [-{hy2m)Vl + U{ri) - fi] 6{ri - r[), 



(D-2) 
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G{ri,r[;i£i), and U{ri,r[;i£i), respectively, and the symbol Tr is defined by Eq. 
(|C-4p without the integration over ri. We then express G{i£i) in the Lehmann 
representatioiPDi''22J as 

Gir^,r[;ze,) = ^ T ^^''"'"''''^ de, (D-S) 

where A{ri,r[;e) is the spectral function. Substituting Eq. ()D-3p into it, we may 
regard Eq. ()D-2p as a functional of j4(e) with the property: 

'" 0, (D.4, 



6A{r[,ri;£ 

which originates from Eq. (|C-5p . Using the procedure to obtain Eq. ()C-10p from Eq 
(|C-8p . we further transform Eq. (|D-2p into 

n = -P — /(e)Tr<^ Imln[^+ Z:(e_) - e_l] 



+Im^(e_)ReG(e_) + Rer(e_)ImG(e_) \ + ^, (D-5) 

where P denotes Cauchy principal value to remove the pole of the Bose distribution 
function in Eq. ()C-9p at 2; = 0, e± = e it iO^ with 0+ as an infinitesimal positive 
constant, and ReG(e_) and ImG(e_) are defined generally as 

R,0(,_) , S£.)±£(£±). Img(..).g'^-':.g<^^>. (D.6) 

respectively. Note ImG(e_) = \A{£) from Eq. (|D-3p . 

We next express ^ of Eq. ()B-17p with A and /. Let us substitute Eq. (|D-laP 
with Eq. (|D-3P into the first-order term of Eq. ()B-17p and carry out integrations over 
r's and summations over Ps with Eq. ()C-9p . We thereby obtain an expression of the 
Hartree-Fock contribution to ^ as 

x[A{ri,rr,e)A{r[,r[-e')±A{ri,r[-e)A{r[,rr,£')]. (D-7) 

Similarly, the second-order contribution in Eq. ()B-17p is transformed into 

^(2) = _ 1 I dr, I dr[ I dr, J dr', y(n - r[)V{r2 - r',) j^^ ^ j"^ ^ 

/OO J />oo J^ 

^ / -^[A{r,,r2;ei)A{r2,rv,e2)A{r[,r',;e',)A{r'„r[;e',) 
"OO "/ — 00 

±A{ri,r2; £i)A{r2, r[;e2)A{r[,r'^; 4)A(r^, n; 4)] J(ei, e'^, £2, 4)> (D-8) 
where J(ei,e'^,e2ie2) signifies 

P Jo JO 
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+(1 o 2) 

+ (1 o 2) 



(l±/l)(l±/()/2/^ 



£1 + e'l - £2 



„-/9(ei+£',-e2-e'2) 
1 + 



/3(ei + e'l - £2 - £2) 



(1 ± /i)(l ± /()/2/^ - /l/((l ±/2)(l±/^) ^ ^j^g^^ 



£1 + e'l - e2 



'2 



with /i = /(ei), /( = /(e'l), etc. The final expression of Eq. ()D-9aP has been 
obtained with 

(1 ± /i)(l ± /0/2/^e-^(^^+^i-^^-^^) = /i/((i ± /2)(l ± /^). 

Note that J(ei,e'i,e2,e'2) is analytic at ei + e'l — e2 — e2 = 0, as is evident from the 
first expression. Hence, it follows that J can be expressed alternatively as 

J(ei,ei,e2,e2J -F— — — ^ -, ^ ^ , , ; -j—, (U-9bj 

ei + Ci — e2 — £2 ei + ei — e2 — £2 

with P denoting the principal value. 

Let us differentiate ^(^^) and ^'^' above in terms of /(e) and compare the 
resultant expressions with Eq. ()B-16p in the same representation. We then observe 
that the following relation holds order by order: 

<5# 1 

;TV^(e)Re^(e_), (D-10) 



Sf{£) 2n 

which apparently originates from Eq. (jB-lSp .^'^ 

With these preliminaries, we now carry out the differentiation of 5 = —dQ/dT. 
Equation (|D-4p tells us that we only need to consider the T dependence of /(e) for 
this purpose. It then follows from Eq. (JD-IOP that the contribution from ^ exactly 
cancels the one from the third term in the curly brackets of Eq. ()D-5p . Hence, we 
obtain 

— 7^Tr{lmln[^+r(e_)-e_l] + Im^(e)ReG^(e)}. 

Next, we express df /dT = —{d/de)k-Q[—f\nf it (1 it /) ln(l it /)] and subsequently 
carry out an integration by parts. We thereby obtain 

9Rer(e_)\ ^^ ^. ,aReG(e_) 



5 = ^B /^ ^a[f{e)] Tr [^(e) (i - ^^^) + 2Imr(e 



de 

D-11) 
where a and / are defined by Eqs. (j4-25p and (|C-9p . respectively, and we have 
removed P from the integration based on the observation that cr[/(e)] only has a 
logarithmic singularity at e = for the relevant Bose distribution function. 

We now consider the case U{r) = in Eq. ()3-2ap . It then follows that every 
matrix such as j4(ri,ri;e) can be expanded in plane waves as 

^(n,r-l;e) = ^^^p(e)e^P-(-^-<)^ 
p 
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with V the volume. Substituting it into Eq. (|D-11|1 and carrying out an integration 
over r, we obtain 



V 



hk^ 



Spde 

{27rhy 



■<T[f{e)]\Ap{e 



1 



dReZ!p{e- 



+ 2ImZ'„(e_ 



dReGp{e^ 



''"^'-' de 

(D-12) 

Noting ReGp(e_) = ReG^(pe) = ReG^(pe) and 2ImZ'p(e_) = As{pe), we reahze 
that this expression is also obtained from Eq. (|4-27p with (j) ^- f. Thus, we have 
confirmed that Eq. (j4-27p is compatible with equilibrium entropy. 
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